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ABSTRACT 


A detailed study of the literature on the concentrations of water 
vapor in the stratosphere and mesosphere is given. It is estimated that 
the concentrations in these regions lie in ibe range from 0. 1 ppm to 10 
ppm. 

A survey was made of the scattering and radiative transfer 
properties of water vapor and the background constituents to determine the 
physical propertier^ of importance to measurements of concentrations. 

It was determined that absorption and emission properties provide 
significant increases in sensitivity compared with the various scattering 
phenomena considered. Microwave absorption in the region of 22 GHz 
and 183 GHz and infrared absorption in the vibrational-rotational band 
systems seem to be the most attractive techniques. Various experimental 
configurations are analyzed and compared. 


vi 



CHAPTER I 


INTRODUCTION 

The problem of how much water vapor exists in the stratosphere 
and above is one that has received attention for a considerable number 
of years but which has still not yet been fully answered. 

The amount of water vapor in the upper atmosphere is a critical 
problem since it plays an important part in the chemistry and heat 
budget of the atmosphere. Water vapor may be responsible fox 
destroying ozone and contributes to the escape of hydrogen through 
photodissociation. It is thought to be responsible for hydrated ions, 

OH airglow emissions, and may control ionization loss processes as 
high as 80 km. Water vapor is involved in O - H - N photochemistry 
in the mesosphere, and there is a strong possibility that the existence 
of noctilucent clouds is dependent on the amount of water vapor present 
near the 80 km level. 

Mastenbrook (1971) reported the observed increase in mixing 
ratios of water vapor from two parts per million to three parts per 
million through the six year period from 1964 to 1969. This increase 
was observed at all pressure levels from 15 km to 21 km. 

For these reasons and more, it is important to develop a method 
of detecting and measuring the concentration of water vapor in the 
atmosphere from the surface of the earth up to about 100 km where the 
H 2 O mixing ratio is probably seriously reduced through dissociation by 
solar ultraviolet radiation. 

At the present time, the highest experimental measurements of 
water vapor have been about 32 km. Thus, there is a large height range, 
30 km to 90 km, where little is known about its concentration. This 
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report discusses various detection techniques that could be used for 
observing water vapor and making measurements. In order to set the 
background, Chapter ' of the report discusses the present knowledge 
about the concentration of water vapor in the atmosphere, both experi- 
mental and theoretical. Also, the attenuation of radiation by water 
vapor is discussed in order to understand the interaction of electro- 
magnetic radiation with the water vapor molecule (which is the basis 
of several of the detection techniques considered). 

In Chapter II tne concentration of the atmospheric constituents 
is discussed, and all that is known currently about the water vapor 
distributions is reviewed. The theory of the various physical radiation 
processes present in the atmosphere is discussed in Chapter III with a 
physical description of the radiative transfer of energy through an 
atmospheric volume, and the refractive index of the atmosphere is 
related to the atmospheric properties of pressure, temperature, 
density of constituents, etc. 

Chapter IV then follows with special attention given to the 
microwave frequency range. The various absorption coefficients for 
water vapor and oxygen molecules valid in this range are presented. 

The theory of infrared absorption is presented in Chapter V 
with the band models of Elsas ser and Goody being discussed. 

In Chapter VI the various methods of radiometry v/hich can be 
used to detect and measure the water vapor concentrations are discussed 
in terms of their basic strengths and weaknesses. The basic configurations 
in which the methods can be used are also discussed. Some of the factors 
which are inherent in the atmosphere and that complicate the various 
detection and measurement processes are considered. 
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In the last chapter an a^terapt is made to draw some conclusions 
about the various methods of detection, the optimum frequency ranges 
of operation, etc. The reasons for the selections are presented although 
at thii; aoint no method seems to be outstanding overall. Each met'^od 
has its weak and strong points. 

The physical and molecular properties of water vapor a 
discussed in Appendix A so that some understanding of the behavior - 
the interaction between the molecule and the incident electromagnetic 
field can be obtained. 

In Appendix B the properties of the oxygen molecule are 
discussed since non- resonant aLsorpLon by oxygen interferes with 
the water vapor lines in the otherwise relatively clean region of the 
microwave spectrum. 

In Appendix C, series approximations to the Mie scattering 
coefficient are presented. 

Finally, in Appendix D, the derivation of Elsas ser's infrared 
absorption model is presented. 



CHAPTER II 


CONSTITUENTS OF THE ATMOSPHERE 


2. 1 General 

The major constituents of the atmosphere up to aboutlZO km are 
molecular nitrogen and molecular oxygen* Also pjresent are small 
amounts of argon, carbon dioxide, water v 2 por, etc* (see Figure 1)* 
Molecular oxygen is dissociated by solar ultraviolet radiation in the 
mesosphere and thermosphere. The concentration of O^, thus, decreases 
while the concentration of atomic oxygen (O) increases with increasing 
altitude. At about 120 km the concentration of O is greater than that of 
at 500 km practically no remains. Molecular nitrogen is much 
more difficult to dissociate and even at 500 km the concentration of 
atomic nitrogen (N) is very small. Above 100 km the effects of 
diffusive separation do become important, and thus, the lighter gases 
tend to predominate in the upper regions compared with the heavier 
gases in the lower regions. As a result of molecular diffusion, heliui 
and hydrogen are the principal constituents in the extreme upper 
atmosphere. 

The concentration of carbon dioxide in the atmosphere is 

approximately constant at least up to 80 km. The maximum o 2 .one 

density is usually found between 20 and 30 km. The number density of 

12 -3 

ozone at its maximum is of the order of 5 x 10 cm . As this means 

5 

there is only one ozone molecule to about 10 oxygen and nitrogen 
molecules, we can expect that ozone will not play a significant part in 
ary nor -re sonant scattering experiments. In the case of water vapor, 
experimenta' information concerning the variation with height above the 
tropopause is poor. 
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2. 2 Water Vapor Content 

The first measurements of stratosj aerie water vapor were 
made by the British in 1953 using a manually operaLad frost point 
hygrometer on an aircraft (Mastenbrook, 1968), The frost point 
hygrometer when used as a sounding instrument provides a vertical 
distribution of water vapor in terms of the frost point temperature. 
From the known relation of frost point temperature and partial pressure 
of water vapor formulated by Goff and Gratch (1946), the concentration 
of crater vapor is determined using the perfect gas law 




or 


Ph,o = '•" 2 ° 


RT 


-H2O 


( 2 . 1 ) 


( 2 . 2 ) 


where q is water vapor density and q is the molecular weight 
2 2 
of water vapor. When Pj^ q, the partial pressure of H^O, and n^j q, 

2 2 
the concentration of H^O, are the saturated values, the value for T is 

the frost point temperature (or dew point temperature where appropriate). 

The British observed in their measurements that humidity 

decreased rapidly with he?ght above the tropopause and mixing ratios^ 

as low as 2 X 10”^ (two parts per million) were found (Mastenbrook 1968). 

Using the same technique, Mastenbrook obtained a median vertical 

distribution of water vapor for the lower stratosphere to a height of 

28 km which approximates a constant mixing ratio within the range of 


^The mixing ratio is the ratio of tlie mass of water vapor present 
to the mass of dry air presen*- in a unit volixme. 
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2x10 ^ to 3 X 10 ^ (two to three parts per million). He found that 
the vertical gradient in the lower stratosphere reversed seasonally. 

The gradient is negative in the late winter and spring and positive in 
the late summer and fall. The change of gradient is attributed to the 
seasonal change of moisture in the low stratosphere (Mastenbrook, 1968). 

A few years ago, Sissenwine (1968) reported that the British 
consider th* mixing ratio to be about two parts per million at the 16 km 
level. In this case the water vapor density tends to become monotonically 
lower as the atmospheric density falls off with altitude. Then since 
frost point (and dew point) have a one to one relationship with absolute 
humidity (water vapor density), they would also have a parallel 
decrease with altitude. Sissenwine (1968) conducted his own investigation 
using a balloon at 32 km and after solving the water vapor outgassing 
problem, concluded that his findings reinforced the British findings 
of a mass mixing ratio of 0. 002 gm/kg at the 16 km level (100 mb) 
which was the maximum altitude of the original British aircraft data. 

His results ranged from 1 to 10 parts per million. These results also 
showed that above this level the water vapor density increased to a 
relative maximum at 25 km (from a humidity minimum above the 
tropopause at about 15 km) (see Figure 2). 

Keneshea, et al. (1972) presented a density profile of water 
vapor which was calculated theoretical'*; ^rom 40 to 150 km and for 
30^N equinox and 60^N wint-^r sing a photochemical model. Sechrist 
(1972) collected theoretical water vapor profiles of Shimazaki and 
Laird (1970), Anderson (1971), Bowman, Thomas, and Geisler ( 1970). 
Shimazaki and Laird and Bowman, Thomas, and Geisler used 
photochemical transport models while Anderson deduced water vapor 
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Figure 2: Water vapor mixing ratio profile from 0 to 32 km altitude 
(Sissenwine, N. et al, 1968; reprinted by permission 
of the American Meteorological Society) 
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concentrations inferred from rocket borne fluorescent scattering 
observations (see Figure 3). 

In one of the more recent papers, Mastenbrook (1971) published 
the results of a six year study in which monthly samples were taken 
over Washington, D. C. and a steady increase of water vapor mixing 
ratio from E x 10*^ to 3 x 10 ^ has been observed. The increase was 
observed at all pressure levels from 100 mb (15 km) to 50 mb (21 km). 
This increase suggests that stratospheric water vapor may possibly 
be an indication of climatic v<iriations occurring at lower levels. 




Figure Theoretical models of water vapor density 
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CHAPTER III 

ATMOSPHERIC RADIATION PROCESSES 


3. 1 General 

The atmosphere is a turbid medium composed of atomic and 
molecular gases and aerosols which are solid or liquid (or a combination 
thereof) particles of different sizes suspended in the air* 

Electromagnetic energy incident upon a volume of this gas- 
particle mixture is both scattered and absorbed in accordance with 
Beer*s law which governs the extinction of energy by the atmosphere. 

where in general 

(T = K + K + K (3.2) 

= the scattering coefficient 

= the absorption coefficient 

K = tlie reflection coefficient 
r 

In addition absorption and emission of radiation occurs in the atmosphere 
entirely independently of external sources of electromagnetic radiation. 

The latter is caused by thermal (black body) radiation resulting from 
the collisions between molecules and particles undergoing thermodynamic 
processes. Thermal or blacK body radiation has a continuous spetrtruim 

ich is described b> Plahck function. Thus the total chknge of intensity 
of radiation in passage through an atmospheric region governed by an 
equation of the form 
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dz 


- cr (I - E) 


( 3 . 3 ) 


where E is the source function related to the thermal radiation intensity 
given by Planck's function 


E oc 


2hv^ 

2 

c 


1 



( 3 . 4 ) 


The physical picture one may draw from this equation is that a 
wave of intensity I incident upon an atmospheric volume in 
thermodynamic equilibri'om is, in general, scattered in directions other 
than that of the direction of propagation (although some forward scatter 
will be in the direction of the wave), reflected by discontinuities in the 
gas mixture, and absorbed. Absorption occurs whenever the electro- 
magnetic energy of the wave raises the molectiles to excited states, 
and then before the molecules can reradiate, collisions occur during 
which non- radiating transitions (deactivation) take place. Emission 
occurs both continuously from thermal processes and discretely 
from molecular transitions triggered by the incident wave. Thus, the 
intensity of radiation leaving the vulume (transmitted) is that portion of 
the wave that has not been scattered, reflected, or absorbed, plus 
that amount which has been forward scattered and emitted by excited 
molecules in the reference direction. 

In order to avoid confusion, the processes should be defined on 
a precise basis. For this study emission is defined as the radiation 
process in which there is a decrease in internal energy of the molecule. 
Actually, there are two types of emission. When the transition to a 
lower state takes place in one step so that the emitted quantum has the 
same frequency as the absorbed quantum, coherent emission occurs, 
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which is also known as resonant scattering. The other type is ca..led 
fluorescent or incoherent emission. In this case the molecule returns 
to the ground state through a cascade process among intermediate 
levels. When frequencies are emitted lower than that of the absorbed 
quantum, these lines are called Stoke'.s lines and the process is known 
as Stoke' s fluorescence. When higher frequencies than that of the 
absorbed quantum are emitted, the lines are called anti-Stoke's lines 
and the process is called anti-Stoke's fluorescence. Anti-Stoke's lines 

_3 

are weaker than Stoke's lines by a factor of 10 because the energy 
in balance must be made up by translational kinetic energy (K.E. ), 

A process related to fluorescent emission is that of Raman 
scattering. Raman scatter differs from fluorescent emission in that 
for fluorescence, the incident photon is completely absorbed «nd the 
molecule is thereby raised to an excited state, whereas in Raman 
scattering, the photon as a whole is never absorbed, but rather simply 
perturbs the molecule and induces it to undergo a vibrational or 
rotational transition and that portion of the photon energy remaining 
becomes r^disttioTi* 

In addition to the process of emission is that of absorption in 
which the energy is transferred to the molecule to raise it to a higher 
state. However, before the molecule can reradiate, the energy may 
be removed in the form of kinetic energy by collisions. All of these 
processes involve a change in the internal energy of the molecule and 
must be described by quantum mechanics in terms of probabilities 
called Einstein coefficients. 

In addition there exist electromagnetic scattering mechanisms 
such as elastic scattering and reflection. Scattering can be defined as 
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an elastic process in which there is no internal energy change of the 
molecules. The process does not emit radiation, but only results from 
a perturbation of the wavefront by the scattering particle. The 
perturbation then acts as a secondary wave propagating out from the 
particle or molecule. This process is described by classical electro- 
magnetic field theory. 

The same is true for the process of reflection. The main 
difference between this process and that of scattering is that scattering 
occurs from individual particles where i >> \ (t being the distance 
between particles and is proportional to n , and X is the wave length 
of the incident wave). However, for the case of reflection, i << \ so 
that a layer of particles creates a fluctuation or discontinuity in 
refractive index, which then causes reflection of the wave. 

3. 2 Emission and Absorption 

Molecules can exist in either the ground state or excited states. 
Though they are neutral, they can possess both electric and magnetic 
'^ipoles or higher order multipoles. The mechanism of energy transfer 
between matter and r. '.Lation is most commonly by coupling between 
the electric component of the radiation field and the electric dipole of 
the molecule or between the molecular magnetic dipole and the magnetic 
field component. This energy transfer may occur over a broad band of 
frequencies (non- resonant absorption) or it may be a sharp line 
resonance transfer (resonant absorption). If, as a result of the 
interaction, the radiation field loses energy to the molecule, absorption 
occurs. Conversely, if the radiation field gains energy from the 
molecule, emission occurs. 
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There are three basic types of radiation field molecule interaction; 
absorption, induced or stimulated emission, and spontaneous emission. 

In absorption a photon is absorbed by the molecule which is then raised 
0 a state of higher energy. Induced emission is the reverse process, 
where the molecule (already in an excited state), under the influence of 
the applied electromagnetic field, emits a photon at the frequency of 
the applied field and drops ^o a lower internal energy level. In 
spontaneous emission, a molecule in an excited energy state, even 
though it is undisturbed by an applied field, emits a photon at the 
characteristic frequency given by the Bohr frequency condition 
(hv^j. = - Ej) and drops to a lower energy state. The processes of 

absorption and induced emission are coherent processes since the same 
incident wave acts on all the molecules, where as the process of 
spontaneous emission is usually an incoherent process since it is 


random in nature. Thus, the resulting spontaneous emission signal 
power is much less than that obtainable in a coherent absorption or 
emission process. The rate of spontaneous emission occurring from 


an upper level E^ to a lower level Ej depends only on the number of 
molecules in the upper level, N^. The rate of stimulated emission and 


absorption depends on the populations and N.. 


These types of 


interaction were dealt with by Einstein. His B coefficient gives the 


probability for either absorption or induced emission, and his A 


coefficient gives the spontaneous emission probability (Barrett, 1958). 


Several life times are important for resonance absorption and 


emission. The excited molecules reach rotational equilibrium before 
they decay back to the ground state because of the fast rotational 

_Q 

relaxation time of the order cf 10 seconds at standard temperature 
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and pressure (STP). The non- radiative vibrational relaxation time is 

-4 -6 

much longer, and is approximately 10 to 10 seconds at atmospheric 
pressure, depending on the molecule. The spontaneous radiation deca> 
time is even longer being about 10”^ to lO"^ or more seconds for 
metastable states. Since the non-radiative vibrational decay time is 
much faster (at STP) than the radiative decay time, the resonance 
emission efficiency is reduced by about the ratio between the two decay 
times (Kildal, et al. , 1971). 

To discuss the process of absorption and induced emission, it is 
helpful to consider a molecular gas having two energy levels Ej and E^, 
having n^ and n^ molecules per cubic centimeter in the upper and lower 
energy levels, respectively. Let be the resonant frequency defined 
by these levels, where hv = E, - E . Then radiation incident on this 
gas in a range of frequencies close to will be absorbed by molecxiles 
initially in the lower state. The incident radiation will supply energy 
hv^ necessary to cause a molecule to make a transition from the lower 
to the upper energy level and a photon of energy hv^ will be removed 
trom the radiation field. Likewise, the radiation can induce or 
stimulate transitions in the reverse direction, from the upper to the 
lower energy level, and this stimulated emission will add a photon of 
energy hv^ to the radiation field. The net absorption by the gas is 
the difference between the number of absorptions and number oi 
emissions. If a molecule in the lower state is capable of absorbing 
energy of frequency v, the probability of absorption per molecule 
per second is B^j p(v>) where p(v) dv is the energy density of the 
incident radiation in the frequency range between v and v -f dv and 
B^j is the Einstein B coefficient for absorption. The number of 
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molecules per cubic ceru. meter capable of absorbing energy at 
frequency v can be written as n^ f(v) where f (v ) is the fraction of 
molecules able to absorb at frequency v. The fraction f ( v) will depend 
on the molecular velocities in the case of a Doppler broadened line, 
on the width of the energy levels in the case of a line broadened h 
stimulated emission, and on pressure in the case of a line broac,ued 
by collisions. In either case, f{v), the line shape function, is subject 
to the condition 




f(v) 


dv = 1 


(3.6) 


o 

Thus the number of absorptions per cubic centimeter per second at 
frequency v is 


n^f( v) Bqj p(u) (3.7) 

Similarly, the number of induced or stimulated emissions per cubic 
centimeter per second is 


nj f(v) Bjq p(v) (3. 8) 

where Bj^ is the Einstein B coefficient for induced emission and is 
related to Bqj by the relation (Barrett, 1953) 

8o®01 = «l®10 

The numbers g^, gj are statistical weights of order unity and are pure 
numbers. They represent the number of sublevels which compose the 
energy levels 0 and 1, respectively (Barrett, 1958). 
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The net number of absorptions is 


"o ®01 ■ 

= n^ £{v) Bqj p(v) 
= Hq f(v) Bqj p(v) 


®10 

- "1 *'oi 

1 °i^° ' 

L ■ “o«i J 


(3. iO) 


Before continuing with the case of absorption and induced emission, 
it is important to discuss the case of spontaneous enroelon a-i ‘■his 
process has been neglected in the net absorption ecuation 


As has been mentioned, the line width can be broadened by 
spontaneous emission or natural broadening. This is because one or 
both of the energy (states) of the transition may have a short lifetime. 
Molecules in such a level or state will not remain in that state if they 
can "make a transition by photon emission". If a molecule is 
assumed to be in a particular state at time t = 0, the probability of 
finding the molecule in that state at time t is the Einstein A coefficient 
and is given by quantum mechanics as 


A = 


64 I M, I ^ 

3hc ^ 


(3. 11) 


where p, is the electric dipole moment. By Fourier analysis of 


radiation from an oscillator whose amplitude is decaying writh time 

-At -At 

according to e energy is decaying at the rate e ), one ci 

obtain a frequency spectnun of the form (Barrett, 1958) 


(v - Vq) + (A/4tt) 


I(v) 


(3.12) 
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or the radiatio^i has a line width A/Ztt. The othe^ forms of line 
broadening, pressure broadening and Doppler broadening, will be 
discussed later in this section. 

Note that A depends on the frequency of the incid, vave v in 
a cubic fashion. Thus, A is larger for optical frequence '.han it is 
for radio frequencies. 

The mechanism of spontaneous emission is the principal 
cause of the continuous emission spectrum generated by the thermal 
processes within the atmospheric gas. The net number of emissions 
is given by 

Iijf(v) A + Hj f(v)BjQ p(\^ = njf(v) [a + Bjg (L,i/)] (3. 13) 

The ret number of absorptions is given by 


nof(v)Boi p(v) 


(3. 14) 


At (thermal) equilibrium, the sum of the emissions is eqvial to that of 
the absorptions 


njf(v) [a + Bjq p(v)] = n^f(v) Bqj p(v) 


,3. 15) 


Then using 


One may obtain 


n. 


n^ 


hv 
o 

1 " kT 

— e 

o 


P(»') iBq 


g, -hi/ g, - hi' 

B -L e Ae — °- 

^10 ® kT J gQ kT 


(3.16) 


(3, 17, 


and it can be shown (Barrett, 1958) that 
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8o®01 


g,B 


I'^IO 


( 3 . 9 ) 


so that 


-hv 

^ ■ ^Tt" 


] p(v) = — e 


-hv 

o 

kT 


or the radiation density p(v») is 


p(v) = 


gi/ggA 


hv 

Tt 


1 ) 


but 


and if gj = gQ 


A = yi h B 

^ 3 01 


p(v) = 


81 ^ V- h 


hv 

o 1 
kT ■ ^ 


( 3 . 18 ) 
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( 3 . 21 ) 


This is the Planck ITack body formula which describes the continuous 
spectrum from diermal raoiation. Returning to the discussion of 
absorption and induct d emission, the relative number of molecules in 
levels 1 and 0 can be related by the Boltzmar. distribution for a medium 
in local thermodynamic equilibrium (collisions dominate radiative 
transfer in determining populations). 


^1 - !l ' 8l -'“•o 
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( 3 . 22 ) 


This can be written for radio frequencies up through the microwave 
range where hv^ << kTas 
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The change in energy density dp(v)dv per unit path length df through the 

h V 3 

gas is — times the net absorptions per cm per second or (Barrett, 1958) 


n»f(v) B-, p(v) h^v^ 


>(w)dv _ -Q i 

di ckT 


(3.24) 


The absorption coefficient K (v) (the fractional change in energy 

a 

density per unit length) becomes 


dp(v)dv 

■ p(v)di " CkT 




(3.25) 


The Einstein B coefficient to be inserted into the above 
formulation is, from quantum mechanics. 


.3, ,2 


(3.26) 


In most atmospheric attenuation work, the total absorption Is 
given in terms of opacity or optical depth t(v) where 


T(v) = J K^(v)di 


(3.27) 


It can be shown that for N molecules in level 0 contained in 

o 

2 

a column along the line of sight of one cm cross section where 




(3.28) 


8ir^ Nq IpI^ 
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(3.29) 


For a line broadened only by spontaneous emissiop, the line 
shape £{v) is (Barrett, 1958) 


r, * 1 Av 
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(3.30) 
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which is proportional to Eqviation (3. 12). Using this expression for 
the opacity of the gas, one can describe how the presence of the 
molecular gas will affec the intensity of radiation received by a 
radiometer (Barrett, 1958). 

Although the preceding theory considered only two states, 
more than two possible energy states exist for all molecules. However, 
it is not possible to cause transitions between every pair of internal 
energy states; restrictions are imposed on what pairs of states may 
be coupled by both radiative and non- radiative transitions. These 
restrictions are known as selection rules and are expressed in terms 
of allowed changes in the quantum numbers specifying the internal 
states. In the ti.eory of interaction of an atom or molecule with electro- 
magnetic radiation^ the effectiveness of radiation on the atom or 
molectile is given by the dipole matrix element p. ^ which is obtained 
quantum mechanically, 'fhis matrix element depends on the initial 
and final states of the transition and like its classical counterpart, 
can be classified as either electric or magnetic. As can be seen from 
either Einstein's A or B coefficients, when a transition takes place, the 
intensity of emitted radiation is proportional ^o the square of the 
matrix element and this provides me key to the selection rules. When 
the matrix elements are calculated quantum mechanically, many cases 
turn out to be zero and it becomes possible to predict quite generally 
which combinations of initial and final state quantum numbers will give 
results different from zero. 


One important property of these transition*’ is that the electric 
dipcle moment is of the order of 100 times the magnetic dipole moment, 
and since the intensity cf emitted radiation varies as the square of the 
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dipole matrix element, electric dipole transitions are some 10 times 
as intense as magnetic dipole transitions. Thus, the concentration 
of molecules possessing electric dipole transitions can be much 
smaller than those of magnetic dipoles and still be detec ^ible. 

There are three main sources of the broadening of the atmospheric 
spectral lines* Pressure or collisional broadening exists because of 
collisions and longer range interactions between the radiating molecule 
and other molecules. In this form of broadening, the atoms or m^’ecules 
are assumed to have a definite collision radius and each radiating 
molecule periodically suffers a sudden disruption of the radiation 
mechanism. The collision either stops the radiating process or changes 
the phase by an arbitrary amount. Fourier integral analysis of the 
wave train yields the observed frequency distribution of spectral 
intensity. The result is tiie dispersion line form (Foley, 1946) 

I(o>) = 3 5* (3-31) 

'' (mr a v) + (a> - 

2 — 

with a half width in angular frequency units of 2mr a v where n is the 
number of molecules per cubic centimeter, a is ttie collision radius 
and V the mean velocity of the molecules. Note that this expression 
is similar to Equations (3. it) and (3. 3U). 

This form of broadening is dominant in planetary or solar 
atmospheres. It is dominant in our atmosphere up to at least the 
upper stratosphere. 

The next most important form of line broadening is Doppler 
broadening which becomes dominant where pressure broadening 
recedes in importance. This process resiilts from the fact that an 
atom or molecule emitting or absorbing radiation is not stationary, but is 
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constantly in motion. For the case of molecules moving along the 
line of sight with velocity v, the freCj lency of its transition will appear 
shifted by an amount given by the well known Doppler formula (Barrett, 
1958) 


V - 



(3.32) 


where v is the observed frequency and is the resonant frequency of 
the stationary molecule. The velocity is assumed to be given by the 
Mavwell'.an velocity distribution 


f(v) 


4tt 


2 1 



■ 2kT 
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where M is the mass of the molecule. For this case, 
distribution of emitted radiation is given by (Barrett, 


(3.33) 


the intensity 
1958) 


I( v) = exp 




which yields the line width as 

V 


= 




tnZ 


(3.34) 


(3.35) 


where is the full width at half maximum intensity. The last 
important form of broadening is that of natural broadening by spontaneous 
emission which has already been described. 

Thus far the absorption and emission processes have been 
general as to the frequency range. The internal energy of a molecule 
is normally characterized by three types of energy; 1) electronic, 
associated with the motion of electrons in various molecular orbits, 
with energy changes typical of visible radiation, *2) vibrational, due to 
nuclear vibrations, with energy changes typical of infrared radiation, 
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and 3) rotational, associated with the end over end rotation of the ent’re 
molecule, with energy changes usually typical of microwave radiation. 
Thus, the absorption and emission spectrum can be conveniently broken 
into various frequency ranges: microwave, infrared, visible, and 
ultraviolet. The most common ranges of interest in studying water 
vapor are the microwave and infrared ranges, although visible 
frequencies are, with the advent of high power lasers, becoming 
increasingly important. For this study only the microwave and infrared 
frequency ranges are considered. 

A problem also to be considered is that in the microwave range, 
pressure broadening differs from that at infrared and optical frequencies 
because the broadening results mainly from non -adiabatic collisions 
(collisions which produce sufficient kinetic energy to make a transition 
from the ground state to an excited state as well as causing a change 
of phase of the oscillator). At infrared and optical frequencies there is 
insufficient kinetic energy available to make this transition and the 
collisions aie adiabatic, the broadening being produced by a change in 
phase of the oscillator due to inter -molecular interactions during the 
collision, rather than a change in energy of the oscillator (Townes 
and Schawlow, 1955). 

3. 3 Raman Emission 

T^e Raman effect is the scattering of electromagnetic radiation 
by matter with a change of frequency. The process is similar to 
Rayleigh scattering where no change in frequency occurs. The 
mechanism of Raman emission can be explained using a classical 
approach, although there are some detail** which must be obtained 
using quantum mechanics. 
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For an incident electromagnetic wave of frequency v^, tlie 
induced dipole moment in a polar molecule may be written as 
(Anderson, lv*l) 

M = ckK cos 2 ttv t (3.36) 

^ o o 

where 

ct is the molecular polarizability 

is the amplitude of the electric field. 

The polarizability will vary with the configuration of the molecule, 
that is, be a function of the vibrational or other motion of the nuclei, 
so that for a given normal mode of frequency 

ft= a +a, co8[2ttv t + ^ (3.37) 

where P is an arbitrary phase. 

Thus, the dipole moment can be written as 

u = E a cos 2 tt V t + E a cos 2ir v t [ cos (Zir v t + p ) ] (3. 38) 

^oo c ol o m 

The molecule can now be viewed as an oscillating dipole which 
radiates pow'er at a rate given by 

1 4 4 _ Z 

I ^ ^ - ||jl| (Anderson, 1971) (3.39) 

3c 

In terms of the preceding eqtiation the expression becomes 

I = E; cos^ 2ir v^t + aj cos^ [ 2tr(v^ + v^)t + ^ 

3c 

2 2 2 

+ o E cos [ 2ir(v - v )t - P] + cross terms 
1 o ' o m' ‘ 


(3.40) 
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Thus, it can be seen, since the cross terms can usually be 
neglected, that the intensity of scattered Radiation contains a term of 
the same frequency as the incident wave and is of the same phase. 

Thus, there is coherent Rayleigh scattering present at all times. The 
other two terms show that there will be a sum and a difference 
frequency radiated with different amounts of phase. The phase will be 
different for each scattering molecule so that the observed radiation 
which is the sum of that from all the scatters will be incoherent. 

This form of scattering provides a unique molecular signature. 
However, it is the weakest of ill atmospheric radiation processes 
considered, having cross section near that ol Rayleigh scatter. 

As wab shown above, for an incident wave of frequency v^, 

Raman lines occur at a series of frequencies v ± v, , v ± v^, ... , 

The lower frequencies are called Stokes lines and the higher frequencies 

are called anti- Stokes lines. The anti- Stokes lines are normally the 

weaker of the two. The frequency displacement is related to the 

rotation- vibration spectrum of a molecule, howeve., selection rules 

differ from those operating in the infrared spectrum (where AJ □ ±1). 

The Raman cross sections vary widely for the individual atmospheric 

components, however, they are in general less than Rayleigh cross 

sections. Thus, to obtain similar signal levels when using Raman 

scatter, some of the range capability available when using Rayleigh 

scatter is traded off to compensate for the smaller interaction 

parameters. The vibrational (vibrational- rotational) cross section 

- 3 

is presently estimated to be about 2.0 x 10 of the Rayleigh scatter 


cross section and a similar situation exists for 
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As was mentioned, the classical analysis does not yield all the 
answers about Raman scatter* Quantum theory explains why the 
anti-Stokes line is weaker* This is due to the fact that anti-Stokes 
lines arise only from excited molecules which are in a minority 
according to the Boltzmann distribution* 

The intensity of Raman lines after averaging over all orientations 
of the molecule and where the molecule makes a transition from state 
m to state n (Derr, et al. , 1970) is given by the following, 


I 

mn 


7 ? ^ 

I 

4 o o 



r>. 

where p, cr are independently x, y, or z. 
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(3.41) 


and 
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(3.42) 


(3.43) 


c = speed of light 

= incident plane polarized intensity of frequency 
r = an intermediate state 

= transition moments of the dipole moment operator 

P- 


There is also the process of resonant Raman scatter in which 
by choosing the incident radiation frequency close to an allowed 
electronic transition of the molecule, an enhancement in the 
scattering cross section occurs due to a resonance effect in the 
polarizability tensor. This enhancement is approximately (Kildal, 
et al. , 1971) 
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(3.44) 


where 

ojj = the pump frequency (incident frequency) 

= the electronic transition frequency 
^ = the line width of the electronic transition 


3.4 Rayleigh and Mie Scattering 

In scattering, the characteristics of the individual scattering 
elements and the statistical properties of the turbid medium greatly 
influence the intensity and geometric nature of the scattering. The 
geometrical structure of the. medium is mainly characterized by two 
dimensionless parameters a/\ and i/\ where a is the average 
diameter of the particle, i is the distance between particles, and 
X. is the wavelength of the incident radiation. The laws of scattering 
differ considerably for cases when a << X and when the particle size 
is comparable to or larger than the wavelength (a^ X). 

When < >> X, the particles can be considered as independent 
radiators. The particle positions are assumed to be random, thus, there 
will be no amplitude interference between the fields of individual 
particles, but only svimmation of scattered radiation intensities. 

The most general scattering theory is that of Mie in which the 
particles are considered to radiate as higher order multipoles 
(quadrupoles, etc.). However, in the approximation when they can be 
treated as radiating dipoles, the theory of Rayleigh holds quite well. 

The classical electromagnetic theory of scattering was first 
developed by Lord Rayleigh and applies to the case where: 



-30- 


1, The dimensions of the scattering particles are much 
smaller than the wavelength (a << X), 

2, The refractive index of the particles differs only slightly 
from that of the medium, 

3, The particles scatter radiation independently of one another, 
this is only possible for i >> X. 


Rayleigh scattering from the molecular atmosphere is important 
for it provides a method by which atmospheric densities may be derived 
from radar or lidar measurements* In addition it also provides a 
convenient, predictable backgroxind to which other scattering and 
absorption effects may be related. For wavelengths well separated from 
the absorption lines of the atmospheric constituents, the Rayleigh 
scattering cross section of an individual scattering center is given 

by (Collis, 1968)* 


ray 


8ir I 
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2ir| 

X 


6 + 36 
6-76 


(3.45) 


where 


6 = depolarization factor due to the anisotropy of the 
atmosphere 

a = molecular polarizability of scatterer 


For the atmospheric gases, the factor 6 has a value near 
0.035, therefore the fraction 

= 1*061 (3.46) 


Note the angular dependence has been integrated out. 
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The polai izability is approximately 2x10 m and thus 

C = 3.96 X 10"'*'^ cm^ (3.47) 

ray 

where \ is expressed in centimeters for the microwave region (X = 1 cm) 

C = 3.96 X lO"^*^ cm^ (3.48) 

ray 

while at the ruby wavelength (X = 0.694 jx) 

C = 1.71 X 10”^^ cm^ (3.49) 

ray 

The total scattering coefficient of a pure gaseous atmosphere 
is the scattering cross section per molecule mult. plied by 

N, the number density of molecular scatterers. 

Thus 

= N C (3.50) 

This quantity is also called the Rayleigh attenuation 

coefficient. When is multiplied by the incident power density 

and the effective illuminated volume, the total power scattered in all 
directions from the incident radiation is obtained. For pure Rayleigh 
scattering, it can be shown that 3/8 ir per steradian of this total will 
be scattered back toward the source. Thus, the volume backs catte ring 
coefficient for Rayleigh scatter is 
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ray 

180 " 8it 


j^ray 

s 


(3.51) 


Consider the following illustrative example. Let a transmitter 
be radiating power P,p from an antenna of gain G (power gain over 
isotropic). This will produce a power density (watts per square meter) 
at a distance Rj (meters) of 
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(3. 52) 


Most bodies, particles, and molecules subtend a small enougn *.ngle 

at the transmitter to be considered point targets. If the point target 

scatters isotropically with an effective cross section C (square meters), 

s 

the scattered power density (flux) at a ra-^ge R 2 from the particle is 


P,pG 

4tR? 


tt S 


47T 


(3.53) 


\ ^2 

The power (watts) returned to the receiver, intercepted by 
an antenna of effective area A (square meters) is (ignoring attenuation) 


PtGCsA 

IbTT^RjR^ 


oc S A 


(3.54) 


For the case of an electrorragnetic wave with electric field intensity 
incident upon a small dielectric sphere of refractive index n, the time 
averaged radiant flux of Rayl^’gh scattering can be written as 
(Kondratyev, 1969) 


I6R^N\^ 


(3. 55) 


where 


E ”P r* 

o oc T watts 


2q 


>1 2 
4ttR m 


power radiated 
2 


m 


S = the scattered energy per unit time per unit area perpendicular 
to the direction of scattered radiation at a distance R from 
the scattering particle at an angle 0 with the incident 
radiation nnpolarized 
c = the speed of light 


q = the (wave) impedance 
of the medium 
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The Mie theory is a comprehensive treatment of the scatteri g 

of a plane electromagnetic wave by a dielectric sphere. The expressions 

are quite complicated involving associated Legendre polynomials, 

Riccati-Bessel functions, and spherical Bessel functions, Howe^'er, 

many useful approximations exist (Appendix C). Mm theory must be 

used whenever the particulate matter has dimensions of magnitude 

similar to the wavelength of the incident ladiation. For large 

particles the elementary scattering cross section C ^ is of the order 

of twice the geometrical cross sectio:.. The scattering pattern in the 

Mie case does not resemble che symmetrical dipole pattern of 

Rayleigh scatter, but can be quite irregular and complicated due to the 

higher order multipole terms. In addition to the index of refraction, 

the results depend critically on the shape and size of the particle. 

Even a sphere has a complex angular dependence. In the case of a 

-3 

water droplet with a radius of 10 mm, the 3cacter is strongest in the 
forward direction, decrea.sing by a factor ^ " t 9 = tt/^ an'^ 
increasing slightly in the direction of 0 ~ 7hus, this case, 
forward scatter and direct ^ ackscatter are the most favorable observing 
angles, Mie scattering from a single particle doet, preserve the initiaJ 
frequency and phase of the incident radiation. However, in an 

ensemble of particles which are (normally) randomly distributed, there 
is a loss of coherence. 

The polarization of the scattered radiation is a complex function 
cf the size parameter = 27ra/\ where a is the particle radius and X 
the incident wavelength, but the degree of polarization is generally 
smaller than in Rayleigh scatter. 
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The ratio of the backscattered intensity to the total scattered 
intensity is a highly variable function of the particle size to wavelength 
ratio and the dielectric characteristics of the particle. Generally, 

Mie scattering is predominantly forward so that in an assemblage of 
particles of different sizes, K' is often less than one in the relation 


B 


Mie 

180 


K' K 


Mie 


(3.56) 


where an approxima^on to is piesented in Appendix C. 

s 


Since the effects of particle cize differences tend to average out 
in reasonable volumes, useful approximate values can be determined 
for K' and used in evaluating the return signal. Calculations performed 
at Stanford Research Institute (Collis, 1968) for water droplet distributions 
typical of natural water clouds give an average value of K = 0. 625. 

Thus, frcm scattering theory one can detect the presence of cloud 
and haze layers, measure their height, shape, and in the absence of 
e 'ressive attenuation, their thickness. Figure 4 shows the difference 


be ween and a clear standaid atmosphere. The magnitude 


of Mie scattering generally depends on the aerosol concentration and, 
thus, varies with changing atmospheric conditions in contrast to molecular 
Rayleigh scattering which remains fairly constant. 

The intensity of the backscatter for both types is given by 

R 


Pr 


87tR^ 


exp 


•|-2 y (T{r)drj- 


(3.57) 


Pr is the received power 
P,p is the transmitted power 


where 
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Figure 4: Volume backscatter coefficients for a clear 

standard atmosphere (Collis, R. T. H. , 1968) 
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c is the speed of light 
T is the pulse duratie ' 

R is the range 

Bi8o is the volximo backscattering coefficient at range R 
A is the effective receiver aperture 
(r(r) is the attenuation coefficient 

3.5 Refractive Index 

The bulk physical properties of the medium and of the 
scattering particles are characterized in general by a compj ^ 
refractive index 

m = n - ix (3. 58) 

where n is the real refractive index and x is directly proportional 
to the absorption coefficient. 

The physical meaning of this complex refractive index is 
as follows: A plane wave of angular frequency 

u, ^ Zvt = ^ (3. 59) 

propagating along the z axis, must have the form 

U = Uq exp[i(wt - Kz)] (3.60) 

in order to satisfy Maxwell's equations for time harmonic fields, 
where U is any of the field variables. The wave number K is then 
written as K = K m where K = ~ is the free space wave 

r\ n n i ^ 


number. 
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One may write 


U = exp[i{^t - K^(n - ix)z}] 

= ®xp[~ exp[i (u) t - K^nz)1 


(3.61) 


From this it can be seen that if \ is not zero, there is a net attenuation 
of the wave. Then since radiation intensity I is p^'oport onal to the 
square of the field strength 


1=1^ K^Z) (3.62) 

where , the coefficient of absorption. 

This complex refractive index means that the atmosphere is a 
loss/ medium (absorption taking place) and can be described in terms 
of a complex dielectric constant 


€ = € j « i€^ (3. 63) 

which will give rise to absorption which can be either non-resonant 
or resonant. In the case of non- resonant absorption, the real part 
of the dielectric constant decreases monotonically as the frequency is 
increased from zero to infinity. The classical theory of non- resonant 
absorption is that of Debye. 

Debye (1929) considered the effect of an impressed electric 
field upon the dielectric constant of both non-polar and polar molecules; 
a polar molecule being one with a strong permanent dipole moment. 

He derived the result that polarizability is composed of two effects: 

The first due to the distortion of all die molecules by the impressed 
field and the second arising from an orientation effect exerted upon 
polar molecules (Bean, 1962). 
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Thus, the polarization P(u>) of a polar fluid under the influence 
of a high frequency radio wave is given by (Bean, 1962) 


, M 4tN 

^ _2 ° 

~ €+ 2 p ~ 3 


“o ^ 3kT 


1 1 

1 + iu> T J 


(3.64) 


where 

€ is the complex dielectric constant 
is the molecular weight 
p is the density o-^ the fluid 
is avogadro's number 

is the average polarizability of the molecxiles in 
the liquid assuming no interaction between molecviles 
p is the permanent dipole moment 
T is the absolute temperature 
k is the Boltzmann's constant 

is the relaxation time required for external field 
induced orientations of the molecules to return to 
random distribution after the field is removed 
(jj is Z'JTv where v is the frequency of the external field 


From Debye's analysis one can conclude that for external 
fields with frequencies less than 100 GHz, diat o>t << 1 and thus 
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(3.65) 


For non-polar gases p = 0 and 
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P(w) = 
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(3.66) 
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For gases the dielectric constant is found to be very close to one so 
that € +2=3. Thus 


P(u>) 


- 1 


+ 2 
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(3.67) 


and 


< - 1 = 4irN a 

o o M 


From the perfect gas law 
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one may write 


« - 1 - 4tN^ “o R T T 
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Likewise for polar gases 


< - I = K^f (A +5, 


(3.68) 


(3.69) 


(3.70) 


(3.71) 


where Kj, K^, A, B are constants. 

For a mixture of gases Dalton's law of partial pressures is 
assumed to hold with the result that one can sum the effects of polar 
and non-polar gases and thus 


j t 

i=l g=l 


For the troposphere, one need only consider the effects of dry air 
(non-polar gases) and water vapor (a polar gas) such that 
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(3.73) 
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•where 

is the pressure of dry air 
is the pressure of water vapor. 

Since the real refractive index is defined as 

n = 

and since the above Equation (3.74) is real. Then 

P, P p. 

- 1 = (A + f ) 

Then one can write 

« , - 1 

n = = V 1 + (€ j - 1) = 2 1 

or n - 1 = 

2 

which is true for « j - 1 <v i or (n - 1 << 1). Thus 

K, P, K, P „ 

. 1 d. 2 

“■^“2 T^2 <p<A+^) 

Since in the atmosphere n is very close to one. it is common to 
a more sensitive parameter, N. 

N, the refractivity, is defined as 

N = (n - 1) X 10^ 

and one then has the general expression for refractivity of the 
atmosphere (Bean, 1962) 

^d 5 Pw 

N = 77.6 T— + 72 -^ + 3.75 X 10 


(3.74) 

(3.75) 

(3.76) 

(3.77) 

(3.78) 
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(3.80) 
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This is an accurate expression since the influence of other 
gases, mainly carbon dioxide, cause less than a 0. 02 percent error. 
Assuming normal atmospheric parameters, the above equation is 
considered valid for frequencies up to 30 GHz. The equation relates 
the bulk atmospheric parameters to the microscopic properties of 
the atmosphere which for large numbers of particles make up the 
macroscopic dielectric properties of the atmosphere. 

It can also be written in tne form (Lusignan, et al, , 1969) 

N = 223 p + 334 p + 1.74xl0^p /T (3.81) 


where the density of dry air and \\ater vapor density are 

3 

expressed in K^/vr. and the temperature T in degrees Kelvin. 

3. 6 Reflection 

This process is a wave interaction with the fluctuations of 
refractive index. For a discontinuity in the refractive index of air, the 
coherent scatter which arises can be approximated with a reflection 
coefficient. This process generally occurs when f << \. This form of 
attenuation arises mainly from aerosols since in general, they have a 
higher index of refraction than air, making a sharper discontinuity. 

This reflection mechanism may be found in a region of temperature 
inversion where an enhanced aerosol concentration may exist. The 
layer of air containing the aerosols will then serve as a coherent 
scattering (reflecting) mechanism of incident electromagnetic energy 
due to the discontinuity (Hajovsky, et al. , 1966). This process is 
dependent on the bulk atmospheric parameters of density, temperature, 
and humidity which determine the refractive index, and is much less 
frequency dependent than the processes of absorption and emission. 
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The theoretical signal level of a radar echo received due to 
reflection from a layer of infinite extent is given by Hajovsky and 
Lagrone (1966) as 


P 


R 


647t^R^ 


(3.82) 


where 

Pr = signal level of radar echo 
P^ = power radiated 

G = gain of the antenna above an isotropic radiator 
\ = wavelength of radiated energy 

= power reflection coefficient of layer 
R = range 


The power reflection coefficient P is given by (Hajovsky, 
etal., 1966) 
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_ Az 8ir _ 


-12 


(3.83) 


for a plane layer with refractive index gradient, AN/Az, and vertical 
incidence. 

The expression for a linear refracHve index gradient is 


where 


AN 

Az 


1.5 


°i - M 
n\*zj 
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p 


(3.84) 


C = (mass) density of the aerosols air mixture in 

3 

micrograms per cubic centimeter (pg/cm ) 
p = density of the suspended particle material in 

3 

grams per cubic centimeter (g/cm ) 
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n 


1 


refractive index of the aerosols (Hajovsky and 
Lagrone, 1966) 


This formulation for holds only for aerosols of diameter 
much less than the wavelength of incident radiation, where the 
refractive index of air is unity, and no interaction between suspended 
particles occur. 

The reflection coefficient may also be expressed from Fresnel’s 
formulas. For normal incidence the fraction of incident flux which is 
reflected in a given direction from the interface between two media is 
defined as (Jenkins and White, 1957) 
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m - 1 
m + 1 


2 



(n + 1)^ + 


(3.85) 


where m is the complex refractive index (m = n - ix). 



CHAPTER IV 


MICROWAVE ABSORPTION 


4. 1 General 

The absorption of tiiicrowave energy from a radio wave by 
atmospheric gases is due to the transition from one molecular rotation 
energy level to another caused by the electromagnetic wave. For this 
process to occur with radiation of a particular wavelength, the molecule 
or atom must have an appreciable absorption cross section at that 
wavelength, and the process must have a high probability compared to 
the other processes that are also energetically possible. 

There are two critical parameters involved in each energy level 
transition. The first is the frequency associated with the transition 
energy. This may be determined from a knowledge of the energy 
levels or from direct measurement. However, the energy differences 
are very small for absorption lines in the microwave spectrum. The 
error involved in taking the difference of two nearly equal values is 
large. Thus, for this reason direct measurement of resonant 
frequency is preferable when possible (Straiton, et al. , I960). 

The second critical number is the line width constant associated 
with the frequency spread of the absorption line. This value should be 
determined by direct measurement for the most reliable results 
(Straiton, et al. , I960). Uncertainty as to the value of Av (the line 
width constant) has been a major stumbling block in the theoretical 
calculation of absorption. 

Some of the atmospheric gases giving rise to (microwave and 
infrared) emission and absorption arc: 02 » ^ 2 ^* ^2' ^2^' 

CH^, and NO 2 (see Figure 5), however, the effects of almost all of 
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theti gases are normally negligible due to their low concentvation. 
Molecular absorption in the centimeter and millimeter wavelength 
bands is primarily due to absorption by oxygen and water vapor (see 
Figures 6 and 7). Water vapor and oxygen molecules have permanent 
electric and mag. etic dipole moments respectively, which when 
excited by an electromagnetic wave, oscillate and rotate with many 
degrees of freedom, each associated with a quantized energy level 
hv. Thus, the molecules absorb discrete amounts of energy from the 
wave and are raised to higher energy levels. In returning to a lower 
level, they reradiate energy isotropically, and therefore, the net 
result is an attenuation of the incident wave. This process is called 
resonant absorption. 

Water vapor has resonances at frequencies of about 22 GHz 
(\ “ 1.35 mm) and 183 GHz (\ = 1.64 mm) (see Figures 6 and V). Oxygen 
has 25 closely spaced resonances occurring in the region of 60 GHz 
(\ »0. 5 cm). However, because of pressure broadening, the wings 
of the individual lines overlap to form a continuous region of absorption 
at elevations near sea level (Altshulter, et al. , 1968). 

Oxygen absorption dominates over that of water vapor except 
in the immediate wavelength region of the water vapor line, even 
though the absorption coefficients of the two gases are of the same 
order of magnitude. This occurs as a consequence of the large density 
ratio of oxygen to water vapor iix the atmosphere. 

The shapes and intensities of the water vapor resonances are 
dependent upon atmospheric temperature, pressure, and the partial 
pressure of water vapor, whereas the shapes of the oxygen resonance 
lines are a function only of atmospheric temperature and pressure since 
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Figure 6; Total attenuation one way transmission through the 
atmosphere (Meyer, J. W., i966); reprinted by 
permission of the Institute of Electrical and 
Electronic Engineers Inc.) 
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Figure *i: Attenuation per kilometer for horizontal propagation 
(Meyer, J, W. , 1966; repr>u #*c by permission of the 
Institute of Electrica.' and Electronic Engineers, Inc.) 
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the concentration of oxygen is known (Altshulter, et al. , 1968). 

Absorption (at 15 - 35 GHz) on clear summer days is caused principally 
by water vapor, while in winter days depending on the latitude, the 
water vapor content is so low that oxygen absorption becomes predominant 
(Altshulter, et al. , 1968), Half of the absorption by water vapor and 
oxygen from the surface up to the itratosphere takes place in the first 
3.2 km which comprises 32 percent of the atmospheric mass. Thi'i is 
due to the pressure broadening of the complex of oxygen lines near 
5 mm. Less than one percent of the absorption is due to the atmosphere 
above 17 km (Shimabukuro, 1966). 

In the portions of the microwave spectrum far ’•emoved from a 
water vapor absorption line, the attenuation is primarily controlled by 
the wings of the submiliimeter (far infrared)lines of and H 2 O. The 
effect on an absorption line can be seen in Fig"re 8. In the frequency 
range 50 to 130 GHz, the first 12 lines (line numbers in the order of 
increasing frequency) make a small contribution to the total absorption 
and the line numbers 21 through 76 make approximately the same 
contribution as the first 20 lines. This set of curves is based on the 
assumption that the line widtl'. constant is 0. 1 cm~^ (3 GHz) for all the 
lines (see Figure 9). It is now believed that the line width will vary 
from line to line, thus the curves most likely Will increase more rapidly 
with frequency than is shown. This is so since Benedict and Kaplan 
(1964) point out that there is considerable theoretical and experimental 
evidence that the line width for the various water vapor lines is 
different depending on the rotational states involved. With the imposing 
array of lines vhich influence millimeter absorption, the general effect 
would be that of having a greater increase in absorption with frequency as 
the shorter millimeter wavelengths are approached. 
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Figure 8: Calculated water vapor absorption for , 5 g/m 
with Av/c = 0.1 (Straiton, A. W. et al, 1960; 
reprinted by permission the Institute of Electrical 
and Electronic Engineers, Inc.) 
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Molecular absorption has been analyzed by many theorists. 

The form presented by Van Vleck (1947) using the Van Vleck-Weisskopf 
(1945) equation has accounted most closely for observed absorption 
data. Quantitative prediction the lines has been less successful. 

Most attempts to adapt the Van Vleck-Weisskopf equation to observed 
absorption measurements use a single line broadening constant 
throughout the absorption spectrum of oxygen and another for water 
vapor in spite of experimental and theoretical evidence to the contrary 
because of che scarcity of information about the exact parameters to 
us The Van Vleck-Weisskopf equation considers the broadening 
effect of collisions between molecules. Pressure or collision 
broadening is the dominant process determining the width of spectral 
lines in the range of pressures that occur from the surface of the 
earth to the upper stratosphere. Under these conditionSt the absorption 
calculated at any frequency that does not coincide with a strong absorption 
maximum is influenced appreciably by the wings of many lines and can 
vary widely depending on the broadening assumed for each of them. 

The major assumption made by Van Vleck and Weisskopf was 
in the collision process where they ass\xmed that after the instantaneous 
impact, the oscillator variables are distributed according to a 
Boltzmann distribution. The Van Vleck formulas givt; the absorption 
per unit distance given the pressure, temperature, and water vapor 
content of a sample of the atmosphere. However, uncertainties still 
exist in the absorption line widths which are temperature and pressure 
dependent, and thus, limit the accuracy of the computed absorption. 

At pressures fo\md in the lower atmosphere, the Lorentz (pressure) 
line widths are of the same order of magnitude as their separation. 
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The general quantum mechanical formula for the absorption 
coefficient for microwaves of frequency s; is (Van Vleck, 1947) 


#1 V wN i 1 

= 10 (log^ge) ^ 


3 

r 1 
3hc 


W '"^1 


I 


^-/kT 

J 


(4.1) 


where |j._ is the matrix element of the dipole moment connecting two 
discrete states i, j which have energy E^, E^ respectively. 


v^j is the frequency of the corresponding spectral line 
N is the number of molecules per cubi<: centimeter 
V is the frequency of the incident radiation 


The factor 10 (logj^e) converts into units of decibels per kilometer, 
and the factor v) is the form factor which determines the shape 

of the absorption line and is given by 




V Av Av I 

ij (v. . - v>) 4 Ai/ (v.. -r v) 4 Ai/ I 

J L' ij ' ' ij J 


(4.2) 


Av is the line width and is proportional to the number of collisions per 
unit time. In order to calculate Av theoretically, the collision cross 
section must be known. 

F rom the fact that 


hv. . = E. - E. = -hv.. 
ij I J Ji 


we have 


(4.3) 


V) 




(4.4) 
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The absorption due to a single line in the millimeter or submillimeter 
regions may then be written as 



(V. 



(10 logjQc) 


Nv 

G 


• [exp(- E./kl) - exp (- E./W)] |p..|^ f(w, ....) 


(4.5) 


where 



G 

h, k, c 


absorption at frequency v due to an absorpt.on line 

centered at v. . in dB/km 
^ -E./kT 

y e ^ , the rotational partition function 


Planck* s constant, Boltzmann's constant, speed of 
light 


Then the factor 


exp(- E^/kT) - exp(“ Ej/kD 

may be expanded as 

exp(- E./kT) [ 1 - exp(- [E^ - E.]/kI)] 
= exp(- E^/ki) [ 1 - exp(h(-Wj.)/kD] 
Then for microwave frequencies 
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hv.j « kT 


exp(- hvj./kl) 



(4.6) 
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or 
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Thus 


! \ I^r\^ ^ » StT Nv 

(V, V..) = (10 log,oe) -q 


h V.. 

■rf 


exp(- E./kT) In. p f(v, v^.) 


The factor vv^j f(v, w.j) may be defined as 


F(w, V..) - vv. . f( V, V. .) 

' ij 


Then the modified Van Vleck-Weisskopf line shape is given by 


(4.9) 


(4. 10) 


F(v, V. .) = V Av 
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ij xj 
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F(v, v,.) = 
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(4. 12) 


Gross (1955) derived a kinetic theory line shape factor assuming 
that after impact, the oscillators have a Maxwellian distribution of 
velocities, but that their positions are unchanged. This line shape 
formula takes the same form as the Van Vleck-Weisskopf formxila at 
resonance, but seems to provide better agreement with measurements 
in the high frequency wing. 

The Gross modified line shape is 
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(4. 13) 



-56- 


where 


V?./ Av 
ij 


V = V. . 

IJ 


V = V. . 
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F(v, V..) 
' ij 


V V.. 4Av/[ V + 4v (Ay) ] v» v. . 

ij 
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(4. 14) 


A third line shape proposed by Schulze -Tolbert (1963) is 


F(v, V..) = 
ij 


where 
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4. 2 Water Vapor Absorption 

Barrett and Chung (19^2) refined the expression of Van Vleck 
on the basis of the experimental findings of Becker and Autler (1946) 
and tne theoretical calculations of Benedict and Kaplan (1959). Barrett 
and Chung jive the pressure and temperature dependence of the water 
vapor line vidths as 

= ‘^‘'o < ioif.2$> + .0046p) (4. 17) 


where x was calculated by Benedict and Kaplan (1964) for H 2 O - H^O 
collisions to be 0. 614. Calculations by Benedict and Kaplan (1959) for 
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for the case of H 2 O - N 2 broadening resulted in = 0 . 626 for the 

2 

temperature range of 220 to 300 K. These calculations indicate a value 

of = 0.09019 cm'^ (2.704 GHz) at T = 300°K. Benedict and Kaplan 

(1964) also found that is approximately equal to X . These 

U2 IN2 

considerations indicate then, that the temperature and pressure 
dependence of the water vapor line widths in air may be represented as 
(Westwater, 1966 ) 
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(4. 18) 


N. 


where P is the partial pressure of the mth molecular species and 


m 


B represents the relative effectiveness of the mth constituent over 
m 


H 2 O - H 2 O collisions in producing water vapor broadening. The 

reference temperature is given by T^ and P = P^j q + Pq + Pj^ is 

2 2 2 

approximately the total pressure of the air. 

The expression used by Chung and Barrett for K , the absorption 

cl 


coefficient at a frequency = 22.235 GHz is 
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where 
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j is the energy involved in the transition (J : 5 - 1 
6-5) and has the energy equivalent value of 446,39 


cm ^ (Dennison, 1940) and 


>>^ 1 , 35 ' 


i2 , 


s the square of the dipole moment determinant for 

-37 

the above transition and has a value of 5. 6 x 10 
(King, etal. . 1947). 


esu 


The first term is a single resonant term representing the water 
vapor line at 1, 35 cm (22. 235 GHz), the second term represents the 
absorption in the wings of all lines below 1.35 cm. In the second term, 
the fact that y « for all lines ether than the 1.35 cm line has been 
used to write 


f(y^j, y) 2y(A/)/Tv^^ 


Then using Van Vlecks (1947) expression for G 


G = 0.034 T 


3/2 


(4. 20) 


(4.21) 


they arrived at 


K^(l. 35 cm) = 1. 05 X 10”^® exp(“^^) 


Ay 


Hv - v^) + 


+ V ) + Av -* 


(v + v^) 


, ,«-52NvAy -1 

52 X 10 — JTT 


(4.22) 


It should be pointed out that the equation that appears in Barrett 
and Chung's paper has two misprints; the above expression has been 


corrected. 
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The line width constant was obtained from Becker and Autler's 
(1946) data by Barrett and Chung as 


where 


Av 


2 . 62 x 10 ^ 

(T/318)* 


H 


z 


P is the total pressure in mm Hg 

3 

p is the density of H 2 O in g/m 


(4.23) 


and T is normalized to 318^K since Becker and Autler’s measurements 
were m'^de at this temperature. 

One problem inherent in this formulation is that the nonresonant 
term was derived under the assumption that the half-width Ay was the 
same for all lines. However, it is now recognized that the absorption at 
frequencies far removed from the resonance line exceeds the theoretical 
predictions. Thus, Barrett and Chung increased the nonresonant term 
by a factor of 5 to match Becker and Autler's measurements. 

Croom (1965) used the Barrett and Chung coefficient for an 
indepth study of wate ’^apor absorption at 22. 235 GHz i.nd in a similar 
fashion obtained an absorption coefficient for the 1. 64 mm line. This 
equation is: 


K ( 1. 64 mm) = 6, 46 x 10 
a' 
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where Ay is the same as Equation (4^23). 


(4.24) 

-1 

cm 


Foa Iho above expression, the first term is due to the resonant 
effect and the second term represents an approximation to the effects 
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of the wings of the 13,5 mm line and the numerous sub -millimeter 
H 2 O lines in the vicinity of 1. 64 mm. 

4. 3 Oxygen Absorption 

The absorption of microwaves is •'•Iso caused by 25 oxygen 
lines at about 60 GHz (5 mm) and a single line at 118.75 GHz (2 .'>3 
Under sea level conditions, pressure broadening of the lines near 
5 mm results in their merging together to form one line, however, 
observations from about 30 km should enable them to be resolved. 


oxygen absorption coefficient obtained from Van Vleck's formula is 


K = 

a 


7 , 


u 


-E 


3(2n + 1) e n 


7 

n odd 


(4. 


2 -E 


[vj = {2l(v +)'(r(v ^;v)|h.^+ 2|(v _)(f(v jv)) |hi^_ + F( Je-^ 
n n n n n 

(4. 

where 


= TV. 


r , Av 1 

22 ^ 2 2 1 

ij (v.. + v) + Av 


Lim 


F(v) = V.. 0 [ V. f(v.., v] = 


(V.. 4 v)‘ 
2vAt 


11 ‘ Ij ■' ij' / 2 . 2 > 

•* Tr(v ^ Av ) 


2 _ 43 ^ n(2n 4 3) 


n 4 1 


(4. 


(4. 


and 


2 ^ 4g‘ (n t 1) (Zn - 1) 

^ - n 

n 

2 83^ (n^ 4 n 4 1) (2n 4 1) 

‘’no ° n(n 4 1) 

he 


P = 


4'n‘mc 


(4. 


(4. 


mm). 

The 

^0) 

26) 

27) 

28) 

29) 

30) 

31) 



-61- 


The selection rvdes for the fine structure transitions permit 

transitions of the form (J“n)~^(J-nt 1) and (J * n) — ^ (J - n 1), 

The corresponding frequencies for these transitions are v ^ and v 

n n 

Meeks and Lilley (1963) modified the expression to the following: 


K^(v, P, T> = C,^ 2, 


n odd 


(4.32) 


where 
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The e. ponent in the Boltzmann factor is 
E 

- 2,06844 (4,33) 

and C, = 2. 6742 for K in decibels/km. 

X a 

Meeks and LiUey (1963) formed an empirical expression for Av 
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where =1.95 (MKz/sec)(mm of Hg)~^, ~ 0*75, and P is in units 

of mm of Hg. 

This empirical expression for the collision line width parameter 
can also be written as 

Av = g(h) (|-) (^) GHz (4.35) 

o o 

where is the standard pressure at the surface of the earth (760 mm 
Hg), is a reference temperature (300*^K) and K' indicates the 
temperature dependence. The altitude factor g(h) is 


g(h) 


/ 

gj 0 1 

(h - h.) 

( gl ^ <g2 - gi> h^ - h^ hj<h < h^ 


V 


g2 


h^ < h 


(4.36) 


as can be seen, Meexs a;id Lilley (1963) used values of gj = 0.60 GHz, 
= 1. 17 GHz, hj = ? km, h.^ = 25 km, and K' = 0. 85. 

At alti*-udef. above 50 km, Doppler broadening of the absorption 
lines becomes comparable to ccllisional broadening and must be 
considered. The expression from Townes and Schawlow (1955) for 
the Doppler line width is 


Avjj = 1. 15 X 10'^ t/(T/T^)°‘® GHz 


(4.37) 


The shape of the Doppler broadened line is different from the 
collision broadened line. However for simplicity, it can be assumed 
that both shapes are defined by the equation for the collision broadened 
line and the resulting line width at lea^ *p to 60 km is Av + 

(Meeks and Lilley, 1963). 
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For oxygen absorption, the nonresonant absorption coefficient, 

0. 0064 dB/km, n:\ust be added to the resonant term to obtain the total 
attenuation coefficient. This term accounts for almost all of the oxygen 
absorption at 10 GHz. 



-64- 


CHAPTER V 

INFRARED ABSORPTION 


5. 1 General 

The principal constituents of the atmosphere are nearly 
transparent in the visible and infrared, but the minor atmosphere 
constituents CO 2 . H 2 O, and O^ have strong absorption bands in the 
infrared, and therefore, according to Kirchoff's law, also emit infrared 
radiation from those bands. 

In the infrared region CO^ bands exist at 1. 4 |j., 1. 6 u, 2. 0 )x, 

2.7 p, 4.3 p, 4.8 p, 5. 2 p, and 15.0 p. There are seven principal 
regions of water vapor absorption in the near infrared beyond 0. 9 p. 
These are due to rotation- vibration bands near 0.94, 1. 1, 1.38, 1.87, 
2.7, 3.2, and 6.3 p. Although each of these regions is sometimes 
referred to as a band, more than one rotation vibration band may be 
present in each region. For example, the 2.7 p band of water vapor 
actually consists of the close lying fundamentals Vj and (Howard, 
et al. , 1956). Other bands are due to methane at 7. 65 pm and N 2 O 
at 7.78 pm. 

The absorption (or emission) in the far infrared (50-1000 pm) 
occurs because of rotational transitions of the molecules while the 
higher energy absorption bands in the infrared region of 1 to 20 pm 
occur because the molecules are vibrating, as well as rotating. 

In a gaseous medium a typical infrared absorption band is 
composed of numerous closely spaced narrow absorption lines due to 
the combined vibrational and rotational energy levels of the absorbing 
molecules. This leads to the absorption coefficient changing rapidly 
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with wavelength, and thus, the averaged transmittance over a wavelength 
interval containing several individual lines may not vary ei.pcnentially 
with the quantity of absorber. However, it is possible to obtain a 
tranfimi ssion fnnrHon which relates the average transmittance over a 
small wavelength interval to an averaged absorption coefficient and the 
quantity of absorber. These transmission functions have been derived 
by Elsasser (1938) for an array of lines whose intensity and spacing 
varies slowly with wavelength and by Goody (1952) for an array of lines 
with statistically random spacing and a uniform or exponential intensity 
distribution. Carbon dioxide absorption has a line structure which can 
be represented by the Elsasser model while water vapor and ozone are 
best represented by the Goody model. 

5. 2 The Elsasser Model of an Absorption Band 

Elsasser (1938) considered the case of an idealized absorption 
band consisting of an infinite array of equally intense, equidistant 
lines, each with a Lorentz shape given by (Howard, et al. , 1?56) 


K = ^ 
a T 


Av 


(v - v^) + Av 


(5. 1) 


where the half width, Av, and average intensity, S^, are considered 
constant for each line. Because the pattern is periodi" ’^requency 
(for equidistant lines) the fractional absorption in any perked is the 
same as the fractional absorption for the entire band. If the line 
centers of tl individual lines are located at v = 0, ±d, ±2d, . * . and 
each line has the Lorentz shape, then the absorption coefficient will 
be given by (Howard, et al. , 1956) 

m 

S. 


K 


I ^ 


Av 


n=-® 


( V - nd)^ + 


(5.2) 
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which can be expressed as 


K 


a 


S 

K(x) = 


sinh P 

cosh p - cos X 


(5.3) 


using the Mittag-Leffler theorem where x = and p = 

one assumes ^hat the incident intensity is constant throughout the 
spectrum, thj fractional absorption for the entire band becomes 
(Howard, et al. , 1956) 

A = -^ ^ (1 - exp[ - K(x) W] )dx (5.4) 

-V 


where W is the absorber concentration in atmos cm. The absorption, 
A^, at any frequency in the vicinity of an absorption line can be 
expressed b.~ 

I -K W 

A = 1 - = 1 - e ^ (5. 5) 

V 1 

o 

where is the intensity of the radiation of frequency v transmitted by 

a sample having absorber concentration W, is the intensity of the 

incident radiation of frequency v and is assumed to be constant for all 

frequencies in the vicinity of the absorpt'on line and K is the 

a 

absorption coefficient at frequency v. 

When the lines are far apart compared to their half width 
(d » Av or p is small) one can write 

S p 

K = K(x) = ^ (5.6) 

2d sin (x/2) 


^ m 

An atmos cm is the path length in cm which would contain the 
same number of molecules if the gas were at STP. 
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and integrate to obtain (Appendix D) 

(ttS 

A = erf 


where 


y 

2 ’ 2 
erf(v) = -F- exp(-t )dt 

VTT ^ 

o 


(5.7) 


The fractional absorption of the band can best be represented by 


A 


erf (-j-) 


1/2 


(5.8) 


where 


21TAVS 

f- = (5.9) 

and is a generalized absorption coefficient for the band. 

The restrictions of the idealized band are relaxed and the 
fractional absorption A is considered to be the absorption at the center 
of an interval Vj to containing several lines that are of approximately 
equal intensity. If one assumes that there is no contribution to the 
absorption at the center of this interval from lines outside the interval, 
the expression can be written (Howard, et al. . 1956) 

i(v,. v,)W 

A(vj, 2 ) = erf( - ) (5.10) 


where f(vj, V 2 ) is the generalized absorption coefficient defined 
above and dependf on the average line spacing d(vj, and the average 
line strength v^) in the interval. The applicable quantity 

f(vi, V 2 ) can be determined empirically over the band and is considered 
to be a function of v, varying slowly over the band, the rapid 
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variation of the primary absorption coefficient K from line to line having 

a 

been obliterated (Howard et.al. , 1956). 

. 3 The Goody Model of an Absorption Band 

In 1950, Cowling (1950) took data on the line positions and relative 
intensities for the pure- rotational absorption lines of water vapor from 
Randall, Dennison, Ginsberg, and Weber (1937). He then assumed the 
Lorentz line shape for the individual lines and computed the absorption 
as a function of absorber concentration for six wavelength regions 
(Howard et al. , 1956). 

These absorption curves were so similar however that Cowling 
was led to the conclusion that in atmospheric work complication is 
avoided and small error is involved if a single absorption curve is used 
at all wavelengths (Howard et al. , 1956). 

This led Goody (1952) to conclude from inspection of the diagrams 
of Randall et al. that the only common feature of these ranges is a nearly 
random distribution of line intensities and positions (Howard et al. , 1956). 

Goody then postulated the so-called statistical or random model 
for a disordered band and showed that this model fits quite well the 
earlier computations of Cowling. He considered the fractional trans- 
mission at the center of an interval of frequencies n 6 wide, where n is 
the number of lines in the interval; and 5 is the mean line spacing. The 
interval is considered sufficiently wide so that there is no absorption 
at its center due to lines outside of the interval. It is postulated that 
all arrangements of line positions are equally probable, that i S , the 
line spacings are random and it is further postulated that there is no 
correlation between line positions and line strengths. The line strengths 
are assumed to be given bv the following function (Howard et al. , 1956) 
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P(S 




( 5 . 11 ) 


where P(S )dS is the probability that a given line has a strength in the 
o o 

range to + dS^ and cr is the mean line strength in the interval. The 
individual lines of the band a e assumed to have the Lorentz shape 
(Howard et al. , 1956). 

With these assumptions about the nature of the disordered band, 
one can show that the transmission (at the center of the interval con- 
sidered, when the radiation traverses the concentration. W, of absorbing 
gas along the transmission path (W);(Howard et al, , 1956) is given by 


T(Av, 


W(t) = exp 
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6 ( Av^ + 


WtrAv 



(5. 12) 


Thus the fractional transmission in such an interval of an absorption 
band depends on the mean line spacing, the mean line strength, and the 
half width as well as on the optical thickness of the absorber. This 
random model can be applied to the near ^nfrared absorption bands of 
H^O. The half width Avis assumed to be the same for all of the H^O 
individual absorption lines although it depends on pressure and tempera- 
ture, also the average line spacing is assumed to be the same for all of 
the H^O bands. Then for a fixed value of W, the transmission at a 
given frequency is considered to depend on the mean line strength tr in 
the neighbourhood of this frequency. The mean line strength is thought 
to be a slowly varying function of frequency over the band and can be 
specified in terms of (the value of W necessary to yield a transmission 
of 1/2 at a given frequency). Thus in the Goody model as in the Elsasser 
lodel, the rapid variation over each individual line of the primary 



-70- 


absorption coefficient has been obliterated. The combination of the 
two fiinctions T(Av, Wcr) and (r(v) then describes the transmission process 
for this particular band (Howard et al. , 1956). It has been shown by 
Plass (1958, 1960) that under conditions when absorption at the line 
centers is virtually complete, W and pressure P enter into transmission 
functions only as a product. The transmission function is termed then 
the strong line approximation. When the absorption is small at all 
wavelengths, including the line centers, the transmittance is independent 
of pressure and depends only on the absorber concentration. This is 
the weak line approximation. 

Beyond a band head where there are no line centers, there may 
nevertheless be absorption from the wings of strong lines near the 
middle of the band. Also in spectral regions where lines from more 
than one absorber are present, the combined transmittance is taken as 
the product of the transmittance s calculated for the separate absorbers, 
provided that there is no correlation between the separate arrays of 
lines. 

The weak line approximation is given as (Bradford et al. , 1964) 


T = 



-S W 
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exp — g— 


(5.13) 


where P = 


ZttAv 


and Y = 


S W 
o 

ZttAv 


The strong line approximation is given as (Bradford et ai, , 1964) 
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1) For the Elsasser model 


Z) For the Goody model 


T = 1 - erf (5. 14) 


T = exp -F ^ (5. 15) 


The wing function for the Lorentz line shape is 





(5.16) 



CHAPTER VI 


METHODS OF DETECTION 


6^1 General 

The techniques of microwave spectroscopy and radio astronomy 
offer a ootentially valuable method for studying water vapor in the upper 
atmosphere. It is useful to characterize these techniques which might 
be adopted for water vapor detection as active and passive. The passive 
method involves only a receiver which can monitor the thermal emission 
(spectra) of the water vapor molecules, or measure the attenuation of 
an external source (usually the sun) by atniospheric water vapor. The 
active method is similar to radar, in which a transmitter- receiver 
system is used to excite a region and then monitor the return. In this 
method, the incident radiation may be characterized by amplitude (or 
power), frequency (or wavelength), phase, direction, polarization, and 
coherence. Any of these characteristics may be chang3d in the process 
of absorption, scatter, or reradiation. Thus the receiver can be designed 
to detect one or more changes of the incident radiation and from these 
changes, inferences are drawn about the content of the region. 

There are three basic types of observation in the passive sense: 
measurement of the spectrum of absorption of the sun's radiation by an 
atmosp* eric constituent, remote (with respect to the earth's surface) 
measurement of the emission spectrum of a constituent, and measure- 
ment monochromatically of the zenith angle (or slant range) dependence 
of the absorption or emission of the atmosphere. 

The emission spectrum measurement utilizes the assumption 
that the atmosphere may be considered to approximate a medium in 
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local thermodynamic equilibrium (LTE) at least up to about 70 km. 

AbovC this, departures from LTE occur because the collisional deactivt 
tion lifetime of an excited molecule becomes comparable with its radiative 
lifetime. 

Thermal emission, when observed from any region of the atmos- 
phere is a function of the distribution of temperature and concentration 
of the emitting gas along the atmospheric path under observation. Thus 
an estimate of concentration can be made by observing the sky temperature 
or apparent emission intensity of an atmospheric constituent. In this 
type of measurement, one requires a reference temperature measure- 
ment. 

In the active method, microwave oscillators or lasers are set 
at the molecular resonances of the crnstituteikt to be measured. Since an 
orders of magnitude increase in reemitted signal may be present, the 
return from the given atmospheric constituent may be isolated. Know- 
ledge of the lirxe shape permits probing in various parts of the line and 
thus differential absorption can provide returns that can also yield 
information on concentration profiles. In ge..eral, since the absorbed 
energy must be at a frequency fixed by the molecular characterist’ 
the source or transmitter should be tunable. 

At the present, the main probing techniques just described can 
be used in vari<.us frequency ranges which include microwave, infrared, 
and visible ranges. There a. 3 certain strengths and weaknesses 
associated with each of these frequency ranges for the detection process. 

Infrared and visible frequencies are still largely a fair weather 
tool as they cannot penetrate significant rain, fog, or cloud cover. For 
this reason, it would appear (on an operational basis) that microwave 
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frequencies are superior for overall detection and actual measurements 
from the ground (or at least within the troposphere). 

Some of the factors which complicate the detection problem are: 
absorption and emission from atoms and nxolecules no‘c of the species 
of interest, spatial and temporal inhomogeneities, scattering by air or 
land, clouds, rain, snow, dust, etc. 

The absorption of radio waves in the lower atmosphere is the 
result of both free molecules and suspended particles such as dust grains 
and water droplets occuring in fog and rain. The effect of absorption 
by liquid water droplets becomes increasingly severe with increasing 
frequency. Clouds and precipitation attenuate electromagnetic waves 
through both absorption and scattering processes. Since cloud droplets 
are generally less than 100 pm in diameter, the Rayleigh approximation 
holds at microwave frequencies. Clouds are the most difficult systems 
to model since they contain water in all three states: gas (uncondensed 
water vapor), liquid (water droplets), and solid (ice crystals). 

The relative humidity in water clouds is for all practical purposes 
100 percent. The amovnt of uncondensed water vapor in clouds approxi- 
mately doubles for each lO^C increase in temperature. Thus the water 
vapor content can vary from ^^4 g/m^ at O^C to ^23 g/m^ at 25^C. 

The attenuation rate for uncondensed water vapor of 0.014 clB/km 

3 

at a wavelength of 1, 5 cm for each g/m of watdr apor content (Dicke, 
1946) <^an lead to several tenths of a decibel per kilometer att nuation 
due to the water apor alone. Thus attenuation from water clouds is 
much larger than that of ice clouds. Moreover the losses for both 
cases, art: due principally to absorption. rather than scattering. Thus, 
the attenuation depends more on total liquid water content than on droplet 
size distribution. 
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In the lower 30 ki.i of the atmosphere, both molecular and aerosol 
scatter are observed. Thus scattering from aerosols must be taken into 
account when considering ground-biased active measurements. 

If the diameter of the particlec is much less than the wavelength 
of radiated energy and the echoes from th? individual scatterers are 
assumed to arrive at the antenna at the same time, but with random phase, 
th'=*n the radar echo signal level from a volume of tl.ese indivMual 
scatterers is (Hajovsky and Lagrone, I960 

2 

I NiOf ,6. II 

i 

where 

Ae is the effective area of the antenna 

T is the pulse length 

c is the spewed of light 

th 

is the mmaber of particles per unit volume of i. size group 

th 

is :he average diameter cf the 'aerosols in the i size group 

This formula holds only under the assumptions that: 1) the 

region containing the aerosols is large enough so that the volume 
illuminated by the radar beam is completely filled with the suspended 
particles; 2) the number of particles or aerosols per unit volume i'"* 
the region is considered to be constant; and 3) the particles are cc .» '»rod 
to be perfect spheres. 

The major equipment problem areas include sensitivitVi absolute 
accvracy, spectral response, and directional response 'antenna. 

In addition, each application has cons ,raints of cost, size, weight, power 
requir aments, availability and operati simplicity. 


P^A T C TT 
p =. 
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Receiver sensitivity may be characterized by the receiver noise 
tem^>erature Tj^, the antenna temperature T^, th^ bandwidth B, integration 
time V, and a constant awhich ie usually in the range 1-3. Thus the equi- 
valent rms flvctuations at the receiver input ATrms is (Staelin, 1969) 




AT 




rms 


K 


( 6 . 2 ) 


This equation may be used to estimate an ultimate limit to receiver 
sensitivity, i v.r example, a radiometer in space viewing the earth would 
see an antenna temperature of ~ 300°K. For a noiseless receiver 

g 

(Tp^ = O) with one second integration time, and B = 10 Hz, ATrms = 

0.03 - 0.09 °K depending on the value oi a . If 1*^K sensitivity is 
sufficient, then an intec >*^ion time ar short as 0.004 seconds may be 
used. These limits ha.. ^ alrong bearing on the ultimate spatial resolu- 
tion and coverage which can be obtained by radi<^meijrs in space (Staelin, 
1969: 

Receiver noise temperaures are in the range of several hundred 
degr es at 3 cm wavelength to about 20, 000 at 3nun wavelength. 

At wavelengths longer tlian 3 cm, receiver noise temperatures of less 
than 150°K can be obtained wl'h solid state parametric amplifiers and 
this capability is being extended to wavele''gths near 1cm (Staelin, 1969). 
Absolute acc'iracies of l^K are xeadily obtainable and values of approxi- 
mately 0. 1 - 0.2 can be obtained with great care (Staelin, 1969). 

For amplifiers other than low ncise amplifiers, bandwidths 
ranging from ! to 10^^ Hz can be obtained over most of the microwave 
region «. £ the spectrum while most parametric or other low noise 
amplifiers have instantaneous bandwidths of approximately 200 MHz or less 
(Staelin, 1969). 
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The directional response of antennas varies considerably 
depending upon wavelength and antenna size. The half- power beam- 
width of most antennas is approximately 1. 3^./D where \ is the wavelength 
and D is the antenna diameter. The sidelobe level should be small and 
in general, tho antenna must be custom built. Some applications require 
scanned antennas (This can be done mechanically or electrically (phased 
array) ). The advantages of mechanical scanning include a superior 
multi-fr^ ^uency capability, lower antenna losses, and electronic 
simplicity. Electric il scanned antennas are more compact, need have 
no moving parts, and can scan more rapidly (Staelin, 1969). 

6, 2 Resonant Absorption 

6.2.1 Passive Configuration 

In this method, the integrated water vapor concentration is 
measured along a path through the atmosphere by monitoring the absorp- 
tion of radiation from a natural source (generally the ean). 

The equation descrioing t^is process for a horizontally 
stratified atmosphere and a flat earth is (Meeks and Lilley, 1963) 



where 

T is the antenna temperature related directly to signal and 
b 

noise pov/er observed by an antenna with an infinitesimally 
narrow beam. 

T^ is the brightness temperature of a source external to the 
measurement region (geneially the sun), 
h ic the height of the antenna ^above the earth's surface. 
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<P is the angle between zenith and the direction of observation. 

(t{zz) is the atmospheric attenuation coefficient at slant range x. 

Th^ sensitivity of the absorption measurement improves with the 
amoiint of water vapor in the path up to unit optical depth (dr = 
thus there exists a unique path length for maximum sensitivity. According 
to Croom (1964), this measurement is about 10 times more sensitive to 
small concentrations of water vapor than is the measurement of an 
emission line. This is due to the greater niimber of xinexcited molecules 
to that number of excited molecules. This method is the most direct and 
least expensive solution to the detection problem and can be used to 
obtain water vapor profiles in principle by taking integrated water vapor 
measurements at several altitudes. 

This method and others can be used to make observations from 
the ground, from space, and by line of sight (occiiltation) measurements. 

A discussion of these observations follow. 

6 . 2 . 1 . 1 Groimd Based Observations (resonant absorption- passive) 
Microwave observations are made of crude total water vapor 
content using broad band radiometers and correlating the attenuation in 
certain AX with ground based humidity measurements. The general 
procedure is to observe the signal level on a number of days and to plot 
the signal strength as a function of the measured .vater vapor content of 
the lower atmosphere. The slope of this line provides the water vapor 
losses and the ordinate or Y axis intercept provides the oxygen loss 
(Straiton et al. , I960). Using this data to predict attenuation under 
other conditions si ^fers not only from the difficulties of the spectral 
characteristics of the emitting source having to be assumed or deter- 
mined separatCi., , but also the composition and other properties of the 
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atmcsphere as a fwinc’"ion of altitude must likewise be assumed or 
measured. The separation of the absorption into water vapor and oxygen 
components then depends on these factors whose effects depend in turn 
on the details of the theory used. 

Another experimental method is to measure the total attenuation 
of an external source of radiation by an atmospheric constituent. In the 
lower atmosphere this method has been used to study water vapor, 
oxygen, and ozone although water vapor and oxygen are now more commonly 
studied by emission techniques, which are discussed in the next rection. 
However, in the upper atmosphere where the concentration of water 
vapor is small, the absorption measurements made from a space plat- 
form may be the only method feasible. 

Advantages of both methods are the following: very sensitive to 
concentrations in the troposphere; simplicity; low power consumption; 
easy acessability of data; no need to create radiation. 

Disadvantages of txie attenixation method are the following: frequency 
of receiver and antenna limited to those of the molecular transitions; 
requires accurate knowledge of absorption line strength; lack of 
sensitivity to concentrations in the upper atmosphere due to the bulk of 
earth's atmosphere attenuating signal; requires inversion of an integral 
equation to extract data. Disadvantages of the correlating method are 
that it is a very crude method although simple to c*mploy and lack of 
sensitivity to concentrations in the upper atmosphere. 

6.21.2 Observations From Above the Surface 
(Resonant Ab s or ption- Passive) 

Inis type of observation can monitor the absorption spectra of 
water vapor if the antenna is viewing an external source (sun, moon, 
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star, etc.). If the antenna is not fixed on ^ source, then the ennission 
spectra of the atmosphere will be seen. 

The upward viewing experiment is the most appealing configura- 
tion for measuring water vapor from a rocket or balloon. The radio- 
meter will nr.easure the radiation absorption due to the water vapor along 
the upward slant path. The Intensity of absorption should decrease and 
approach zero as the rocket altitude increases since the concentration of 
water vapor is thought to decrease greatly above 90 km. This decrease 
in intensity when contrasted against the height of the rocket or vehicle 
should determine the relative water vapor concentration. 

The main reason for using an absorption experiment from space 
is due to the enhanced sensitivity to concentrations above 30 km. The 
disadvantages are as follows: the equ'pment must have greater sensitivity 
and frequency stability to monitor the less broadened lines; restrictions 
on volume, weight, and power available; a data transmission system 
must be included; and an elevated instrumentation platform (rocket, 
balloon, satellite, etc.) is required. The main problem for a space 
based absorption experiment is to maintain ' le antenna fixed on the source 
throughout the Aight. Another problem to be considered is that the inte- 
gration time will have to be quite fast in order to use a rocket, yet the 
sensitivity of a radiometer when viewing the solar disc for example, is 
low imless the integration time is long. Thus, there exists a difficult 
compromise between sensitivity and integration speed. 

^ . 2 . 1,3 Satellite Occultation (resonant absorption- passive) 

III occultation, a satellite-borne radiometer receives the radiation 
emitted from another satellite (active resonant ab«-orptiDn), from the 
atmosphere (passive resonant emission), or from the sun (or other 
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natural source) along a ray path (passive resonant absorption). The ray 
path may be identified by the height (Tangent height) of the point (Tangent 
point) closest to the surface of the earth (see figure 10). 

An important quantity measured in the occultation experiment is 
the radio refractivity , This occurs since radio waves propagate in the 
atmosphere not with the speed of light c. but with a reduced speed 
given by 

V = c/n where (6.4) 

P 

n is the refractive index 

The modulus of refraction or refractivity N, is proportional to 
both the air density and to the water vapor density (see equation 3. 81). 

The changes in refractive index have two effects on the radio waves. 

First, the phase velocity of the propagated wave front changes differen- 
tially causing the actual ray path to be deviated from the straight ray to 
a curved ray path (this is due to the atmospheric refractive index gradient). 
Second, the path appears to be longer since the waves travel slower, 
that is, the waves take a longer time to traverse the path and therefore, 
there are more wavelengths along the path. Thus, in its travel over the 
curved path, the wave suffers an overall phase retardation <Ps given by 
(Tank, J969) 


0s 



I 


N(s) ds 


(6.5) 


where 

X. is the radiation wavelength 
ds is an element of the curved path 
N(s) is the modulus of refraction 
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Absorption is another quantity which can be measured in occulta- 
tion since the net effect of the atmosphere on the occuitation signal is due 
to both refractive variations and absorption. In 1966, a Sto ford study 
(Liusignan et al. , 1969) showed that the main source of error in the 
occuitation system is water vapor in the lower atmosphere and that using 
the amplitude of the occuitation signals themselves, the absorption due 
to water vapor could be used to determine the average water content. 

This effect will be superposed over the effect of refraction. Most of the 
absorption would take place within a fraction of a scale height ot che tan- 
gent point since the atmospheric density and pressure fall off exponentially 
with height. The refractive effects will also occur near the tangent point. 

The Stanford group also proposed the use of a passive microwave 
sensor at 22.235 GHz to monitor the absorption due to water vapor and 
thus provide for an independent source of absorption data. To date, the 
use of the passive absorption experiment has had difficulty with vertical 
resolution' (Lusignan et al., 1969). 

'^he occuitation technique became an accepted experimental process 
when the Mariner 4 and 5 spacecraft detected surprising density profiles 
in the martian atmosphere using a single frequency bistatic radar measure- 
ment of the phase pach variations. For more detail on the theory see 
Fjeldbo et al. , (1965). 

Advantages of this measurement technique are the following: 
long term measurerrien* an be made (compared to a rocket or air- 
craft); large areas can oe monitored in a short period of time; the out- 
gassing of water vapor is less of a problem after a relative short time. 
However, equipment packages capable of operating over long periods of 


time must be obtained. 
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6.2.2 Active Configurations (resonant absorption) 

In this method, the integrated water vapor concentration is 
measured along a one way path through the atmosphere usi.ng the absorp 
tion of an active source (a source which requires a separate supply of 
power to create radiation). 

The frequency of the source is first tuned to the center of an 
absorption line. The intensity of radiation at the detector is then given 
by 


I 

on 


I T . 


o air 



K N (R) dR 
a 


( 6 . 6 ) 


whero 

is the intensity of the radiation source 
T^.^ is the atmospheric transmission function for dry air 
R is the path length 

The transmitter is then detuned so that the frequency is outside 
of the absorption line. The intensity received by the detector is then 


off 


I T . 


o air 


(6.7) 


The integrated concentration is determined from the ratio between the 
on and off line intensities and is given by 
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N (R) dR 
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6*2.2 1 Ground Based Observations 
(resonant absorption-active) 

In this configuration a ground based radiometer could measure 
the absorption of radiation from a rocket, balloon, or satellite* However, 
only the total integrated water vapor density will be known. Thus this 
observation will suffer from the lack of vertical resolution. The experi- 
ment has been preformed using a ruby laser emission ^ine at 6943.00 
and tuning the laser cavity in order to shift the line so that it c oincides 
with the 6943.815 line of water vapor. The measurements have only 
been made in the lowest few kilometers of the atmosphere, but the method 
is capable of extension to greater heights. 

The advantages of this system are the following: one has control 
of the source characteristics; noise can be minimized by using narrow 
bandwidth signals. Most all the other ground based equipment advantages 
apply here also. 

The disadvantages are that the frequency of the transmitter must 
now be controlled and that radiation has to be produced. This means 
that permission must be uined from federal regulatory bodies. 

Also, the receiver and transmitter have to be in different locations, 
which can lead to a very expensive experiment. 

6. 2. 2. 2 Observations from Above the Surface 
(resonant absorption- active) 

This configuration could be used between tv;o different space- 
craft. For example, a fixed balloon platform at 30 k.n containing the 
source and then monitor this source with a radiometer cn board a rocket. 
It may also be used between two rockets or two satellites. The latter is 
the acti\e occultatmn experiment and consists ol two satellites in orbit 
which can be spaced so that different Ixne of sight paths can be chosen 
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between the two. This experiment has great promise for examining quite 
carefully upper atmospheric layers. The theory of this type of experiment 
has been discussed in section 6.2. 1.3. 

Advantages of such systems are the ability to control the source 
characteristics and that the measurements would be free of the effects 
of lower atmospheric water vapor. However, the vehicles carry water 
vapor into the measurement region, thus providing a possible error source. 
Another problem to be dealt with is the exist -'e of a shock wave which 
builds up around any supersonically traveling body such as a rocke*-. 

This shockwave is a layer of air at a different temperature than the air 
making up the atmosphere. This difference in temperature gives rise 
to a difference in refractivity between the antenna and the external 
atmosphere. As the speed of the body increases, the temperature of the 
layer increases and at a. sufficiently high temperature, the layer becomes 
a plasma (ionized gas) which can prevent the reception of radio signals. 

However, even more important is the fact that two independent 
vehicles introduce such complexities into the measurement that its 
feasibility becemes vanishingly small. 

6. 3 Resonant Emission 

6.3.1 Passiv Configuration 

This method also the least costly, monitors the emission spectra 
of water vapor from any point in the atmosphere. The equation describing 
this process for a horizontally stratified atmosphere and a flat earth is 
(Meeks and Lilley, 1963). 



(6.9) 
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where 

is the antenna temperature observed by an antenna with 

an infinitesimally narrow beam 

T is the kinetic temperature of the atmosphere 

h is the height of the antenna above the earth’s surface 

<p is the angle between zenith and the direction of observation 

cr(x) is the atmospheric attenuation coefficient at slant range x 

This equation describes the emission taking place within the 

atmosphere along the slant path. 

6. 3. 1.1 Ground Based Observations 
(resonant emission- passive) 

Spectral emission tecluiiques are presently used for the bulk of 
the microwave atmospheric studies. Those observations near a wave- 
length of 5 mm (or 60GHz) are dominated by the thermal emission from 
molecular oxygen. By varying the depth of penetration into the atmosphere, 
the temperature structure in the lower atmosphere may be determined. 

Two methods are available for varying the depth of penetration. First, 
one can make observations at different frequencies at a fixed zenith 
angle, such that those frequencies least attenuated by absorption features 
of oxygen will penetrate the atmosphere further. Second, one can make 
observations at different zenith angles, at a fixtd frequency so that at 
large zenith angles, longer paths through the atmosphere are probed 
than at smaller zenith angles. As the path length thrcugn atmosphere 
increases with zenith angle, there is a corresponding increase in attenua- 
tion for a plane parallel atmosphere. The lath length is proportional to 
the secant of the zenith angle and if the atmosphere were horizontally 


stratified 
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T (<f>) = sec (f} 


16, 10) 


where 

<t> ir the zeaith angle and 

is the vertical atmospheric opacity or optical depth 

For the case of the spherical atmosphere, the secant law in 
ufeful up to zenith angles of about 85° (Wuifsberg, 1967). 

In both of the above methods, it is necessary to have temperature 
information. Also equation 6.9 must be inverted. 

The inversion process is very important to the accuracy of the 
observational data. Considerations for this process are the following: 
how does the technique handle noise in the observational data, h ^ 
fast does the particular technique converge to a solution and is tnat 
solution unique and real. The most serious error is imposed by numerical 
dixferentiation which grossly amplifies any noise or observational error. 

Spe*'tral emission observations at wavelengths near 1.3 3 cm 
detect the thermal emission from water vapor and have been shown to i* ' 
an effective method of inferring crude tropospheric water vapor profiles, 
integrated water vapor, and stratosp'" cric water vapor. The observations 
can also differentiate liquid water from water vapor by virtue of the 
different spectral responses in frhe two cases. For this purpose, a 
microwave channel in the 0. P to 1 cm range of wavelengths where these 
differences are e?sily detectible is desirable. Microwave se ors come 
closer than any other sensor to having all weather capability. Like most 
passive microwave atmospheric studies, the method may be regarded 
as still in the stage of infanc / and its full potential has not been fully 
explc'* d. However, since ground based measure, lents of water vapor 



vre not sensitive to the atmosphere above S kin, it is necessary to 
consider space based observ-itions , 

Advantages of this passive system are that there is no need to 
create radiation and the lew cost of an experiment. Heweve-. this 
experiment is not as sensitive to sr^all concentrations as is that of 
resonant absorption. 

6.3. 1.2 Observations from Above the Surface 
(resonant emission- passive) 

Deteccion of Wevter vapor in the upper <itmospher be 

accomplished using vehicles in space in the following manner, \) 4ew 
the surface of the earth from above in which case the eo rs>vity of the 
surface viewed must be known, this ic discussed further in the text. 

2) view the atmosphere along an upward slant path where the emission 
spectra of the water vapoi present along the line of sight path is moni- 
tored. 

Spaceborne, downward looking experiment;, at wavelen^^tis nea 
the water vapor resonance line at 1.35 cm are capable of sensing atmos- 
pheric water vapor and liquid wate' '.ontent of clouds oceans. This 

method is at least an order of magnitude less sensitive over land sur- 
faces. The major reason for the difference in the two cases is that the 
microwa/e temperature of the oceans is 150^K because of the low 
emissivity of water, whereas the microwave temperat’^re of land is 
~ 300 K. S’noe the atmospheric water ' j,por is near 280 K, the contrast 
between the wate - vapor and h ckground is considerably enhanced over 
the oceans, Staelin (1968) reported that tnis configuration could be used 
to determine the integr-ited wa^er vapor ontent to within abouL 0. ' gm-cir 
when coupled with a microwave temperatui v ounding experiment. How- 
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ever, this temperature soxinding experiment will present additional 
complications in a spa.eborne probe. 

As was mentioned in section 6.2,1. 2, the upward viewing experi- 
ment seems very appealing. For this configuration however the radiometer 
must scan only the old background atmosphere 300^K) in order to 
measure the emission spectra. This experiment is less sensitive than 
the passive absorption measurement. But if sensitivity is not a problem, 
then the emission experiment could be much easier to preform since it 
^hould be T .ach easier to keep the antenna fixed on the cold background 
<han on a hot source of radiation which has a much smaller solid angle. 

Most of the other advrntages and disadvantages which have been discussed 
in section 6. 2. 1 . 2 apply here also. 

6.3.2 Active Configuration 

(resonant emission- active) 

In this method also known as resonant scattering, the molecules 
in a region are excited by the transmitter beam and then the excited 
molecules emit radiation spontaneously into a solid angle of 4ir 
steradians. Measuring the emitted radiation determines the relative 
concentration of the molecules. 

This emission only occurs when the wavelength of the incident 
radiation is close to that of an absorption line. Since the cross section 
will vary across the absorption line and since most experiments 
performed use radiation with a bandwidth greater than that of the absorp- 
tion line, it is convenient to discuss an integrated back- scattering function 


I 

Av 





p (v) dv 
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where P(v) is the backscattering function at frequency w. It can be 
shown that 




e^f ' 


€ me 
o 


( 6 . 11 ) 


where 

^ is the oscillator strength of the transition (the equivalent 
number of harmonic oscillators required to yield the same 
strength) and is given by 

f ' = B ' (6. 12) 

n n 2 r n 

ire 

e is the charge of an electron 
m is the mass of an electron 

is the free space permittivity 
c is the speed of light 

h is Planck's constant 

and 

B ' is Einstein's coefficient for induced emission 
n n 

The value of B(v) will depend on the oscillator strength of the 

transition as well as the line shape, but will be many orders of magnitude 

greater than the Rayleigh scatter. 

6. 3.2.1 Ground Basect Measurements 
(resonant emission-active) 

This method has been put to use in recent years using high power 


lasers to excite the medium. At the present, the state of the art (size 
vs. power) requirements of lasers limit their use to ground based 
measureme; ^:s. This mechanism does have one strength which may 
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make it a very promising technique in the future. Derr et al. (1970) 

reported that resonant erriission is quenched or perturbed and is rarely 

observed at high atmospheric pressures. This is because resonant 

- 8 

emission has characteristic re-emission times of the order of 10 
seconds to 10"^ seconds whereas the average time between air mc^f*- 

_ 9 

cule collisions in the troposphere is usually very short ( ~ 10 

seconds;. Fiiteen kilometers is, with few exceptions, a lower lim'* jn 
the observation of resonant emission (Derr et al. , 1970). 

The method has only been used to detect sodium at very high 
altitudes thus far, however, there is only a small concentration of 
sodium. The detection was made possible by the very strong resonant 
line of sodium. Thus the technique has potential for use in detection 
of any molecule which has a strong resonant line such as water vapor. 

Advantages of this system are that control of the source allows 
accurate phasing of the source pulse to the receiver and narrow band- 
width signals can be used to limit noise. However, a strong power 
source is required. This is no problem on the ground, but in space, 
it is out of the question at the present time. 

6. 4 Raman Scatter 

In this method (always active), a source of fixed frequency is 
used to excite a volume of molecules which in turn will scatter (reradiate 
instantaneously) the energy with a shift in frequency characteristic of 
the molecule. 

The return power from the range R is given by (Kobayasi, 

197v) 


P 

r 


P^ f k T^T^A^Y(R)N(R) o-;. (ir) 



(6.13) 
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where 

is the transmitted power 
i is one half the pulse length 

k is the total efficiency of the receiver and transmitter 

T. is the atmosphere transmittance at the transmitter frequency 

is the atmosphere transmittance at the Raman shifted frequency 

A is the effective receiver area 
r 

Y(R) is the geomefrical factor that accoimts for overlap of the 
transmitter and receiver beam paths 
N(R) is the density of the molecules and 

is the Raman vibrational* rotational backs cattering differen- 
tial cross-section 

The Raman scatter cross-section is inversely proportional to 

the fourth power of the wavelength, thus the largest cross-section would 

be in the ultraviolet region of the spectrum. The visible region has the 

next largest cross-section (~ 10* cm /sr). However, of all the processes 

considered, this me*:hod has the smallest c rons-section (the average 

- 3 

Raman cross-section is of the order of 4 x 10 that of the Rayleigh 

cross-section for atmospheric gases). This is the main problem with 

using this method in the microwave region where the cross-section is 

46 2 

of the order of lO" cm /bt. 

6.4.1 Ground Based Observations 
(Raman scatter- active) 

This method has been used for detection of the major atmospheric 
components such as and O^, also for H^O in the lower atmosphere 
where the density of H 2 O is greatest. This technique has not been used 
(to: the author's knowledge) for observations above the surface (excluding 
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smokestack observations). This is most likely due to the low sensitivity 
of the method to the small concentrations of molecules and atoms present 
in the upper atmosphere. At present the method is limited to ground- 
based measurements, thus space-based observations will not be consider- 
ed. 

Advantages of Raman scatter are in its ability to avoid completely 
the ambiguity in the return at the transmitter frequency due to the Mie 
and Rayleigh scattering components. Moreover, the Raman method 
requires only a single wavelength source where in general, the resonant 
absorption and resonant emission schemes require a tunable source. 

Disadvantages in addition to the small cross-section are the 
following: that since all molecules in the atmosphere contribute to the 
Raman scattering, good filtering techniques are required to select the 
Raman signal for the constituent of interest; the requirement of high 
power transmitters; the lack of sensitivity over long distances (Kildalj 
etal., 1971). 



CHAPTER VII 


SUMMARY AND DISCUSSION 

Present knowledge of water vapor concentrations above 30 km is 
extremely meager „ Estimates of volume mixing ratios of 1-3 ppm are 
widely used, however realistically one cannot rule out uncertainties of 
plus or minus an order of magnitude about this levels These uncertainties 
complicate our understanding of fundamental atmospheric processes and 
structural models. Consequently the significance of an upper atmospheric 
water vapor measurement is becoming clearer. 

The present work is intended to form a comprehensive back- 
ground for the development of a measurement technique. It deals with 
the following areas of interest: what is known and presumed of water 
vapor in the upper atmosphere; the molecular properties of H^O, 
especially in the microwave but to a lesser extent in the infrared; 
radiative transfer and scattering theory as they apply to atmospheric 
measurement problems; and a survey of measurement c .xfigurations . 

In the microwave region, water vapor has two rotational lines 
at about 22 „ 2 GHz (1.35 cm) and 183 GHz (1.64 mm). The 22 GHz line 
is by far the weaker of the two but lies in a relatively uncluttered region 
of the spectrum (see Figs, 11 and 12) for which adequate transmitting 
and receiving equipment is already developed. Millimeter wave 
instrumentation development could easily make the 183 GHz line rather 
attractive for the mesospheric measurement where sensitivity is essential. 
There is a close lying line, as shown in Fig. 12, however it is at 
least an order of magnitude weaker. 
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Figure 11: Zenith attenuation versus frequency {Crane, R. K. , 
1971; reprinted by permission of the Institute of 
Electrical and Electronic Engineers, Inc.) 
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Figure 12: Line frequencies of several atmospheric gases 
(Tolbert, C. W. et al, 1963; reprinted by 
permission of the Institute of Electrical and 
Electronic Engineers, Inc,) 
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There are a multitude of vibrational- rotational bands of water 
vapor through the infrared* Among these, strong bands between l,3p 
and 1. 9|i and at 6. 3 p appear quite attractive and should be looked at in 
more detail. Tuneable laser technology at these wavelengths is now 
available promising rather sensitive absorption measurements from 
active configuration experiments , 

The satellite occultation configuration is clearly superior to the 
others considered. First of all, it alone could provide a truly global 
water vapor measurement which is necessary for the speciiic agronomic 
and ionospheric problems which bear on the stratospheric and mesospheric 
water vapor concentrations. Secondly, the large atmospheric paths 
lengths associated with the occultation configuration appear to be very 
necessary in view of the low water vapor concentrations likely to be 
found in the mesosphere. Finally, maximum sensitivity for microwave 
as well as infrared wavelengths, appears to require an active, two 
vehicle configuration. The simplest configuration which meets all 
these criteria Involves two vehicles in closely spaced station keeping 
orbits. 

The next step appears to be detailed engineering calculations 
with specific hardware and orbital vehicles in mind. 



APPENDIX A 


THE WATER VAPOR MOLECULE 


The water vapor molecule (H^O) is an asymmetric top molecule 
and has a very complicated spectrum. The dipole moment is very strong 
(1,85 Debye) and thus the molecule undergoes very strong electric 
dipole transitions. The molecule has three non- degenerate eigenvibrations, 
two of which are symmetric and one antisymmetric. All three of its 
principal moments of inertia are different. The molecule can be represented 
by a plane model which has the property that the sum of the two principal 
moments lying in the plane is equal to the moment of inertia perpen- 
dicular to the plane. This, however, does not simplify the equations of 
motion. In the unexcited state, the water vapor molecular configuration 
has the form of an isosceles triangle with side length = 0. 958 R 
(nucleus to nucleus) hydrogen oxygen distance and an apex angle 6 = 

104^ 30'. The absorption of radiation by water vapor in the region from 
0. 54 to Op is caused by vibrational- rotational transitions and in the far- 
infrared (from 9p to 1.5 cm) by purely rotational transitions. In the far 
ultra-violet (\<0.2p), the absorption is determined by electronic transi- 
tions, dissociation, and ionization. 

The rotation- vibration spectrum has been accurately mapped from 
135p to 5700 R. Ghosh and Edwards (1956) compiled a listing of H 2 O 
resonance lines in the microwave and far infrared regions. The pure 
rotational spectrum reaches its maximum intensity aroiind 50p, from 
which point it slowly decreases until at about 15p only a few lines are 
observable. At about 8. 5p, an intense absorption again sets in, this 
time due to the o’*ter edge of a fundamental vibration band. 
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The only significant transitions occuring at frequencies up to 
300 GHz are at 22. 235 GHz (J : 5 6 ,) and at 183.31 GHz (J : 2.,^3 _). 

Detailed information on the 22.235 GHz line has been difficult to 
get with conventional absorption spectroscopy because of the s^aU 
absorption coefficient (7.21 x lO"^ cm“^). 

Randall et al. (1937) gives a listing of the rotational lines of water 
vapor which is quite extensive. The selection rules governing fh allowed 
transitions of water vapor are given by Dennison (1931). 

The selection rules follow from the wave functions and their 
symmetry properties which describe the energy l. vels involved in the 
transitions. For the electric dipole, the selection rule is AJ = O, ± 1 
except J = 0-~- = O which is forbidden. For each quantum level J, 

there are 2J + 1 nondegenerate energy levels. These energy levels are 
distinguished by the subscript t where t - - J for the lowest level, 

- J + 1 for the next higher le/el and so on up to +J the highest level. 

A transition from one energy state J^to another .T^ of higher energy is 
accomplished by the absorption of a discrete amount of energy in accor- 
dance with the selection rule. AJ > 1 is forbidden in a dipole transition. 
The transitions are divided into three branches relating to whether the 
change in J (AJ) is -1, O, or +1 and are known as tbj P, O, and R 
branches, respectively. AJ = +1 corresponds to absorption and AJ = -1 
corresponds to emission thus the P branch does not contribute to the 
absorption of energy (Hall, 1967). 

Each level of the molecule is given a symbol ( -*-+), (+-), (-+), 
or (--) which denotes the symmetry group to which the level belongs. 

Since in water vapor, the electric moment lies along the middle axis of 
inertia, the only radiative transitions are those connecting the levels 
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(++) (--) or ( + -) *-*(- +). There are three types of normal vibrations 

of the water vapor molecule which are observed to have the following 
main wave numbers. 

^ O ^ O 


Vj = 3670 cm” ^ = 1675 cm" ^ = 3790 cm’ ^ 

(\j = 2.7p) (X^ = 6p) (\3=2.64^) 

All of these types are active in absorption. 

The rotational water vapor absorption spectrum is associated 
with rotational energy transitions caused by the molecule rotation 
around its axes a, b, c, asymmetrically and in a spinlike manner (see 
figure 13). 

Since the selection rules only allow the transitions (++) to combine 
with (--) and (- +) to combine with ( + -), there is no absorption of the 
non- resonant type since the electric moment matrix of the molecu. 
contains no diagonal elementc. Non-resonant absorption would require 
diagonal matrix elements and thus would demand that (++) combine with 
(++) which is not allowed. 




APPENDIX B 


MOLECULAR OXYGEN 

The oxygen molecule has a number of transitions in the micro- 
wave region around 60 GHz and a transition at 118.75 GHz. The rnicio- 
wave spectrum of the oxygen molecule is a result of fine structure 
transitions in which the magnetic moment assumes various directions 
with respect to the rotational angular momentum of the molecule. 

Oxygen has an intranuclear distance of 1. 128 R. The molecule does 
not possess an electric dipole moment however it does have a magnetic 
dipole moment. This magnetic moment arises because the molecule 
possesses two unpaired electrons. These unpair<»d spins give a spi*' 
quantiim number of S = 1. The rotational quantum number N (N describes 
the end over end rotation of the molecule) can only assume odd values 
because of the exclusion principle. 

An odd rotational state N then combines with a spin state to 
produce a total angular momentum vector with quantxun number J = N - 1, 
N, or N + 1. Selection rules for the fine structure transitions permit 
transitions of the form (J = N) -►(J = N + 1) and (J = N)-*(J = N - 1). 

The gyromagnetic ratio g ^ Lande factor g = 2, Thus the 

oxygen molecule has a magnetic dipole moment of two Bohr magnetons. 
This magnetic moment interacting with the magnetic field of the micro- 
wave radiation results in absorption of energy. At high altitudes, the 
Zeeman effect in oxygen cannot be neglected. The earth's magnetic field 
will split the oxygen lines into separate components. Each quantun. state 
of the molecule specified by the quantu^ii number J will be split into 
2J + 1 magnetic substates. The fine structure transitions AJ = ± 1 
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proceed according to the selection rules AM = O, ± 1 where M is the 
magnetic quantiim number. Generally there are 3(2J' + 1) Zeeman 
components in the transition where J' is the smaller of the two J values 
involved. In the case of oxygen, there is non- resonant absorption since 
the absorption is caused by the magnetic dipole moment for which the 
rule is that + combines with + and - combines with - . This gives rise 
to the matrix for the magnetic dipole moment having diagonal matrix 
elements . 



APPENDIX C 


MIE SCATTER APPROXIMATION 


Penndorf (1962) has obtained an approximation formula which 
will extend the Rayleigh range by about a factor of 2. This formula 
can be derived from Mie's formula using the series expansion method. 
For real indices of refraction n, the total scattering coefficient is 


where 


K 


mie 


8a 


2 ,\2 

n ^1 

n^ + 2 


6a 


1 +■ 


n -2 


n^+2i 


“I 


n^ + 41n^ - 284n^ + 284 
175 1 2 2 

' (n'^ + zr 


(C.l) 




2ira 


(C.2) 


and , the size parameter. Note that the first term is the Rayleigh 

formula 


K = 
s 


128 TT^a'* 


- 1 
n^ + 2] 


(C.3) 


This approximation (C. 1) can be used up to a size parameter of = 1,4 

and n < 2. The error between the exact computation (Mie) and K will 

s 

remain below 2 percent up to n - 1.5 and reach 15 percent at n = 2. 

For the case where the indices of refraction are complex 
(m = n-ix)» the total extinction coefficient is 
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+ 4 . 8n X 




(C.4) 




36 


2 


where 

Zj = (n^ + x^) + 4(n^ - x^) + 4 


Z 2 = 4(n^ + x^) + 12(n^ - X^ + 9 

The total scattering coefficient is 



(C.5) 


(C.6) 

(C.7) 


These two equations (C.4), (C.6) lead to reliable values for ~ 0.8 

in the range n = 1. 25 to n = 1.75 and x ^ 1 • For any larger particle the 
exact Mie forumlas have to be used. The useful size rai - fcr water 
droplets at weather radar frequencies for which con\pl c indices of 
refraction pertain, is restricted to a.>out 0.2. 



APPENDIX D 



through the substitution 


(D.6) 
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Z = sf 2 Z 

'' U - Of 


To evaluate this integral, let 


® 2 2 
r> -ax 


■I* 


w 

I' (a) = -2a y e 


2 2 
-a X 


® 2 2 , 

e* _-a X dx 


dx + 2a 


X + 


= - ,/V + 2a'tt- I (a) 


dl , 2 . - 

Tsr - 2a« I = - <fir 


which has the solution 


2 2 
-a a 
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f «■ 


dt + c 


Then from (D. 8) 


= t Tan 
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which implies that C = , thus 


2 2 2 2 - 2 
T, , •trao' Traof2 l-u 

iw = IT' T^r J ' 


du 


= yI" ^ 1 - erf (a a ) J 


(D. 13) 


where 

erf (aa) 



For the integral in equation (D*6) 
a = 1 and the result is 

A = erf (or) where 

JS WirAv 
ar = ^ o 

d 


(D. 14) 


(D. 15) 
(D.3) 
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